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We  begin  with  an  introduction  to  the  classic  Black-Scholes  option  pricing  model  in 
chapter  1 . This  explains  the  method  by  which  stochastic  processes  are  used  to  obtain 
arbitrage  prices  for  options.  Chapter  2 follows  with  recent  results  in  the  field  and  touches 
on  how  some  of  these  results  will  be  extended  in  chapter  4. 

Chapter  3 provides  some  new  results  on  Brownian  motion,  which  are  of  great  use 
in  chapter  4.  Chapter  4 begins  by  introducing  the  first  new  class  of  options,  the  partial 
barrier  tunnel  options.  It  then  proceeds  to  price  these  options  using  the  results  fi’om 
chapter  3.  Chapter  5 is  a follow-up  of  chapter  4.  It  first  shows  how  the  results  of  chapter 
4 generalize  those  in  the  existing  literature  and  then  provides  numerical  results  and  an 
analysis  that  illustrates  how  changing  the  parameters  affects  the  price. 
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Chapter  6 introduces  another  new  option,  the  Get-Out  option.  This  option  depends 
on  two  underlying  securities  whereas  all  of  the  aforementioned  options  depend  only  on 
one.  The  goal  again  is  to  derive  an  expression  for  the  arbitrage  price  of  this  option.  We 
separate  the  pricing  into  two  cases.  The  first  case  assumes  that  the  two  underlying 
securities  are  independent  and  therefore  do  not  require  joint  distributions  for  the 
processes.  The  second  case,  however,  assumes  that  they  are  correlated.  Consequently,  the 
pricing  is  not  as  straightforward  and  requires  joint  distributions  of  the  processes  and  some 
of  their  functionals. 


vii 


CHAPTER  1 
INTRODUCTION 

1.1  Background 

The  financial  markets  are  growing  more  rapidly  now  than  ever  before.  Investors  at 
all  levels  are  using  financial  products  to  hedge  their  positions,  that  is,  to  help  reduce  their 
risk.  As  a result,  the  derivatives  market  has  grown  as  well.  One  of  the  most  basic  and 
earliest  financial  instruments  used  for  hedging  is  the  standard  European  call  option  on 
stock.  This  call  option  gives  its  owner  the  right,  but  not  the  obligation,  to  purchase  a 
security  at  a prespecified  price  and  time  in  the  future.  The  prespecified  price  is  called  the 
strike  price  and  the  prespecified  time  is  the  expiry.  A put  option  is  similar  in  nature  but 
gives  the  holder  the  right  to  sell  instead  of  buy.  American  options  are  very  similar  to 
European  options  but  they  allow  the  owner  to  exercise  the  option  at  any  point  up  until 
expiry,  whereas  European  options  may  only  be  exercised  at  expiry.  From  here  on  we  will 
work  with  only  European  options. 

Suppose  one  holds  a long  position  in  a particular  stock.  Purchasing  a put  option  on 
the  same  stock  would  guarantee  this  investor  the  ability  to  sell  the  stock  at  the  strike  price. 
The  downside  loss  would  then  be  reduced.  Determining  prices  for  options  similar  to  these 
but  with  additional  stipulations  is  the  point  of  this  work. 
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1 .2  Arbitrage  Pricing 

One  can  see  that  the  price  of  a standard  call  option  must  be  less  than  or  equal  to 
the  current  price  of  the  underlying  stock.  If  this  were  not  the  case,  one  could  sell  the  call 
option  and  buy  the  stock.  Since  the  option  costs  more,  this  investor  would  also  have  some 
cash  left  over.  Furthermore,  at  expiration  only  two  things  can  occur.  If  the  call  option 
were  exercised,  he  would  only  be  left  with  the  cash  left  from  the  original  transaction.  If  it 
is  not  exercised,  then  he  could  simply  sell  the  stock  back  into  the  market  and  have  even 
more  money.  In  particular,  if  the  call  option  is  priced  higher  than  the  stock,  then  there  are 
opportunities  to  make  money  without  any  initial  investment  or  risk.  Thus,  if  one  were  to 
purchase  a call  option,  he  would  not  pay  more  than  the  current  stock  price.  Arguments 
similar  to  this  dictate  ranges  in  which  prices  must  fall. 

Definition  1.2.1  An  arbitrage  opportunity  is  a situation  in  which  one  may  make  a profit 
with  no  initial  investment  and  without  taking  any  risk. 

If  two  exchanges  carry  the  same  security  but  at  different  prices,  then  an  individual  could 
simultaneously  purchase  the  security  at  the  lower  price  and  sell  at  the  higher  price.  This  is 
an  example  of  an  arbitrage  opportunity.  Next  we  will  see  how  arbitrage  can  be  used  to 
price  a call  option. 

1.3  The  Discrete  Case 

1.3.1  The  Model 

We  begin  with  the  simplest  model  available  for  a stock,  the  single  step  binomial 
model.  This  model  assumes  that  a stock  either  goes  up  or  down  by  certain  percentages 
over  the  lifetime  of  the  option  [1].  For  simplicity  we  will  assume  that  we  have  zero  interest 
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rates  here.  Denote  by  Sq,S^,Sj  and  it  the  initial  price,  the  price  if  the  stock  goes  up,  the 
price  if  it  goes  down  and  the  strike  price. 


Figure  1.3:  Discrete  model 

If  the  stock  goes  up,  then  the  value  of  the  call  option  is 

-A:)"  =max(^„ -A:,0).  (1.1) 

If  it  goes  down  then  it  is  worth 

- ky  = max(5^  -k,0).  (1.2) 

One  could  also  choose  to  create  a portfolio  of  the  stock  and  a cash  bond  valued  at  1 dollar 
so  that  the  net  worth  of  the  portfolio  would  be  the  same  as  the  value  of  the  call  option  at 
expiry.  To  see  how  to  do  this  we  start  with  x units  of  stock  and  y units  of  bond. 
Depending  upon  whether  the  stock  goes  up  or  down,  x and  y must  satisfy  the  equations 

xS,  +>;1  -ky 

(13) 

in  order  to  replicate  the  payoff  of  the  call  option.  We  can  solve  this  system  for  x and  y 
so  that  the  portfolio  has  the  same  terminal  value  as  the  call  option.  The  cost  of  composing 
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such  a portfolio  would  then  have  to  be  the  appropriate  price  for  the  call  option.  Any  other 
price  would  lead  to  an  arbitrage  opportunity. 

Derinition  1.3.1  The  arbitrage  price  for  an  option  is  the  value  of  the  portfolio  that 
replicates  the  payoff  of  the  option. 

1.3.2  The  Martingale  Measure 

Notice  that  we  did  not  use  any  probability  in  the  preceding  arbitrage  arguments. 
However,  pricing  the  call  option  in  this  manner  does  indeed  induce  a probability  measure. 
Let  C be  the  price  of  the  call  option  and  Sj  be  the  price  of  the  stock  at  expiry.  We  would 

like  a probability  measure  P under  which  C is  the  expected  value  of  the  possible 
outcomes.  That  is,  P must  satisfy 


(1.4) 


where  /?,  and  p.^.  the  probabilities  of  the  stock  going  up  = 5„)  or  down 

{St  =Sj)  under  P Since  /?,  and  p^  must  also  satisfy  /?,  + = 1 , it  is  clear  that  there  is 

a unique  solution  for  /?,  and  . Moreover,  the  relationship  [Sj  -kf  <C<  {S^  - k)* 


must  hold  in  order  to  prevent  arbitrage  and  this  ensures  that  both  p^  and  p^  are 
nonnegative.  The  probability  measure  P induced  in  this  manner  has  an  additional 
implication.  We  must  first  clarify  a couple  of  concepts.  We  first  introduce  the  more 

general  discrete  process  S = {5’,  defined  on  a sample  space  Q . The  process  S may 

branch  up  or  down,  as  in  the  single  step  model,  at  each  step.  The  sample  space  Q is 
composed  of  all  possible  trajectories  of  the  process.  This  more  general  model  is  known  as 
the  Cox-Ross-Rubinstein  model  [9]. 
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Definition  1.3,2  The  filtration  generated  by  the  discrete  process  S = {5,  is  the 
collection  of  a - algebras  f-a  (5^  : j < ij  generated  by  the  process. 

Definition  1.3.3  A discrete  process  S = {5,  is  a martingale  with  respect  to  a measure 
P and  its  filtration  if  £’|5,|  < oo  for  all  i and  = 5",  for  all  i<  j. 


Under  this  measure,  the  stock  price  turns  out  to  be  a martingale.  In  general,  the  interest 
rates  are  not  equal  to  zero.  When  interest  rates  are  nonzero,  we  must  adjust  for  the  time 
value  of  money.  In  particular,  receiving  x dollars  at  a future  time  t is  equivalent  to 
purchasing  e~''‘x  worth  of  a bond  with  interest  rate  r now.  Note  that  invested  in 
this  bond  will  payoff  e"'  = x dollars  at  time  t . So  the  value  of  receiving  x dollars 

at  a future  time  t is  worth  e~'’'x  dollars  now.  This  is  what  we  call  a discounting  factor.  So 
in  the  case  of  nonzero  interest  rates,  the  discounted  stock  price  process,  e~’^^'S,  , will  be  a 
martingale  under  this  measure. 

Definition  1.3.4  A martingale  measure  is  a measure  under  which  the  discounted  stock 
price  is  a martingale. 

Once  we  have  a martingale  measure,  we  may  price  other  options  by  simply  taking  the 
expected  value  of  their  discounted  payoff  under  this  measure. 

Theorem  1.3.4  The  arbitrage  price  C for  an  option  is  given  by  C = E’^  [X]  where  P is 
the  martingale  measure  and  X is  the  discounted  payoff  of  the  option. 
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We  omit  the  proof  of  this  theorem  as  it  may  be  found  in  [9],  This  reduces  the  problem  of 
pricing  options  to  finding  the  martingale  measure.  Once  the  martingale  measure  is  found, 
we  may  simply  take  the  expected  value  of  the  discounted  payoff  to  obtain  the  arbitrage 
price  of  any  option. 

1.4  The  Continuous  Case 

1 .4. 1 The  Black-Scholes  Model 

We  now  move  to  a better  model  that  allows  for  more  than  just  two  possible 
outcomes  for  the  stock  price.  In  fact,  this  model  allows  the  stock  price  to  take  any  positive 
value.  A stochastic  process  is  used  to  model  the  stock  price  and  an  exponential  function 
(deterministic)  is  used  for  the  bond  price  [1]. 

The  Black-Scholes  model  for  options  pricing  is  the  most  widely  accepted  and  used 
model  in  the  field  [3].  Despite  some  of  its  unrealistic  assumptions  such  as  continuous  time 
trading  and  no  transaction  costs,  the  model  has  persisted  in  both  academia  and  industry 
since  its  birth  in  1973.  We  first  fix  a probability  space  on  which  we  define  a 

standard  one-dimensional  Brownian  motion  W, . Denote  by  T,  = <j {W/,0  < s < t)  the 

a - algebra  generated  by  The  basic  Black-Scholes  model  is  set  forth  by  the  following 
stochastic  differential  equations 

dS,  - oSfdW,  + juS^dt 
dB,^rB,dt 

where  denotes  the  stock  price  process,  B^  the  bond  price  and  cr,  and  r are  constants 
representing  the  volatility,  drift  and  interest  rate,  respectively.  Corresponding  to  each 
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option,  there  is  a payoff  that  is  dependent  upon  what  the  underlying  security  does  by  the 
expiry  time  T . 

Definition  1.4.1  A contingent  claim  X is  an  J-j  measurable  random  variable. 

A contingent  claim  is  used  to  describe  the  payoff  of  an  option.  Thus  the  contingent  claim 
for  a standard  European  call  option  is  X = (5;.  - e . 

1.4,2  The  Martingale  Measure 

Just  as  in  the  discrete  case,  we  need  only  find  a martingale  measure  under  which 
we  can  compute  option  prices  by  simply  taking  the  expected  values.  Though  not  as  simple 
as  solving  a system  of  linear  equations  as  in  the  discrete  case,  we  are  guaranteed  a 
martingale  measure  by  the  well-known  Cameron-Martin-Girsanov  theorem  [7].  Under  this 
measure  the  discounted  stock  price  process  is  a martingale.  In  other  words,  if  we  let 

g 

^ be  the  discounted  stock  price  process,  then  we  have  dZ,  = oZdW.  where  W is 

' ' ' 

a Brownian  motion  under  the  martingale  measure.  Thus  we  may  price  options  by  taking 
their  expected  value  under  this  measure.  Since  this  measure  always  exists  in  our  setting, 
we  will  assume  that  the  measure  P we  use  to  price  options  is  the  martingale  measure. 
Furthermore,  we  will  drop  the  W,  notation  and  use  as  a Brownian  motion  under  the 

martingale  measure  P An  important  consequence  of  using  the  martingale  measure  is  that 
the  original  drift  term  fi  drops  out.  In  order  for  the  discounted  stock  price  process  to  be  a 

martingale,  the  drift  term  in  the  process  must  be  r - — [7].  Since  the  original  drift  term 
is  irrelevant,  we  choose  to  start  with  the  martingale  measure  and  avoid  it  altogether. 
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1.4.3  An  Example 

Take  a standard  European  call  option  on  a stock  with  volatility  cr , drift  ^ . Let 
k and  T be  the  strike  price  and  expiry,  respectively.  Let  P be  the  martingale  measure. 
Under  this  measure  we  have 


fits,  = <rS,dW,  + r S,dt  (1.6) 

where  W,  is  a Brownian  motion  under  P . Note  that  the  drift  //  is  irrelevant  since  we  are 

using  the  martingale  measure.  The  solution  to  this  stochastic  differential  equation  is  given 
by 


S,  = exp 


^ 2 

N 

afV,  + 

a 

r 

t 

V 

^ 2 J 

y 

(1.7) 


We  need  only  compute  C - e (5'j.  - A:)  J . We  use  the  superscript  to  emphasize 


that  this  expected  value  is  taken  under  the  martingale  measure  P . This  expression 
depends  only  on  the  marginal  distribution  of  Sj- . After  making  appropriate  substitutions 


we  may  express  the  price  C in  terms  of  the  normal  cumulative  distribution  function  as 


C = S^N 


In  ^ 

V ^ J 

+ 

( 1 2^ 

r + -a 

1 2 J 

■> 

T 

1 

r 

1 

f 

In 

—2. 

. k J 

+ 

f 1 2^ 

r--a 

1 2 J 

N 

T 

crVr 

V j 

ayff 

^ J 

(1.8) 


where  N (•)  is  the  normal  cumulative  distribution  function  (i.e.  N(x)  = P(X  < x)  where 

A'  is  a standard  normal  random  variable).  Again  notice  that  the  drift  term  n does  not 
appear  in  this  formula. 


CHAPTER  2 
RECENT  RESULTS 


2.1  Barrier  Options 


One  of  the  first  modifications  that  can  be  made  to  the  basic  call  or  put  option  is  the 


see  whether  or  not  the  underlying  stock  hits  the  barrier  during  the  course  of  the  option. 
These  new  options,  called  barrier  options,  come  in  two  types:  knock-out  and  knock-in  [9]. 
A knock-out  barrier  option  would  be  the  same  as  a regular  option,  but  with  the  added 
stipulation  that  the  stock  price  must  not  hit  the  barrier.  If  the  stock  price  hits  the  barrier 
during  the  course  of  the  option,  then  the  option  automatically  expires  with  zero  payoff. 


On  the  other  hand,  a knock-in  option  has  zero  payoff  by  default.  The  option’s  payoff  takes 
place  only  if  the  stock  price  hits  the  barrier  during  the  lifetime  of  the  option.  For  example, 
the  payoff  of  a knock-out  barrier  call  option  could  be  of  the  form 


addition  of  a barrier.  A barrier  is  a prespecified  price  level,  which  is  monitored  in  order  to 


h 


0 


T 


Figure  2. 1 : Barrier  option 


(2.1) 


9 


10 


where  h is  the  barrier  level.  In  this  example  the  barrier  is  placed  below  the  initial  stock 
price,  that  is  h <Sq.  This  type  is  commonly  referred  to  as  a down  and  out  option.  Pricing 
such  an  option  depends  not  only  on  the  marginal  distribution  of  Sj- , but  also  on  whether 
or  not  the  stock  price  hits  the  barrier.  Depending  upon  whether  the  barrier  is  below  or 
above  the  current  stock  price,  the  joint  density  for  the  pair  ) or  {SjMt  ) is 

necessary  [12].  Here  m,  and  M,  denote  the  running  minimum  and  maximum  of  the 
process  W,  as  given  by 


m,  = inf  W. 


0<s<t 


M,  = sup 

0<s<J 


(2.2) 


These  densities  are  readily  found  in  most  texts  on  stochastic  processes  including  [7]  and 
[10].  The  density  is  given  by 


P (fV,  e dx,Ur  ^dy)=  exp 


^ i2y-xf^ 
2t 


dxdy . 


(2.3) 


The  barrier  option  price  admits  a representation  by  a linear  combination  of  cumulative 
normal  density  functions  [9], 

Barrier  options  have  a variety  of  applications  in  the  markets.  Consider  the  case 
where  one  has  a future  cash  flow  in  another  currency  [2].  The  exchange  rate  affects  the 
value  of  this  cash  flow  relative  to  the  investor.  A down  and  out  call  option  on  the 
exchange  rate  would  enable  him  to  hedge  his  position.  If  the  exchange  rate  fell  sharply, 
then  the  call  option  would  be  of  little  use  anymore  and  the  option  get  knocked  out.  Thus 
the  introduction  of  the  barrier  would  make  the  call  option  less  expensive  for  the  investor 


and  still  suit  his  needs. 
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2,2  Partial  Barrier  Options 

A further  generalization  of  the  barrier  option  is  the  partial  barrier  [4],  The  only 
difference  between  a partial  barrier  option  and  a regular  barrier  option  is  that  the  barrier  is 
not  monitored  for  the  entire  lifetime  of  the  option.  The  monitoring  period  would  either 
start  at  time  zero  and  end  at  some  time  before  expiry  or  it  would  begin  at  some  time 
before  expiry  and  end  at  the  expiry.  Notice  that  for  the  latter  case,  the  boundary  could  be 
hit  from  above  or  below. 

A 

h- 
0 

Figure  2.2:  Partial  barrier  option 

These  two  cases  are  distinguished  and  priced  separately  in  [4].  The  payoff  for  a partial 
barrier  option  where  the  monitoring  period  starts  at  time  zero  is 

^{S,^AV/€[0,/,  ]}  (2.4) 

where  <T . Conditional  expectations  and  the  Markov  property  are  used  in  [4]  to 
compute  the  expectation  of  each  payoff.  The  prices  are  expressed  as  linear  combinations 
of  bivariate  and  univariate  normal  cumulative  distribution  functions.  The  bivariate 
cumulative  distribution  function  is  defined  as  N{X, Y, p)  = P{^X  <x,Y  <y)  where  X 
and  Y are  normal  random  variables  with  correlation  coefficient  p . A more  direct 
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approach  for  several  of  the  computations  above  is  found  in  [2],  The  joint  density  for  the 
pair  ) where  5 r is  derived  and  then  used  to  price  a knock-out  call  option  where 

the  partial  barrier  begins  at  time  zero.  Not  surprisingly,  the  joint  density  is  a linear 
combination  of  bivariate  and  univariate  normal  density  functions. 

Partial  barrier  options  may  be  of  use  both  individually  or  when  packaged  with 
other  options.  The  nature  of  this  option  allows  one  to  compensate  for  changes  in  volatility 
over  the  lifetime  of  the  option.  One  may  expect  moderate  volatility  in  the  near  future  but 
increased  volatility  soon  after.  In  this  case,  one  could  choose  a partial  barrier  option  with 
the  barrier  only  monitored  in  the  near  future.  This  would  serve  to  decrease  the  cost  of  the 
option  and  still  serve  to  hedge  or  speculate  an  underlying  security.  One  may  also  package 
barrier  and  partial  barrier  options  together  to  create  ladder  options  [4].  A ladder  option  is 
just  like  a regular  call  or  put  option  except  for  that  it  locks  in  intrinsic  value  when  certain 
barriers  are  reached. 

2.3  Double  Barrier  Options 

Another  modification  made  to  the  standard  barrier  option  is  the  addition  of  a 
second  barrier.  Naturally,  the  two  barriers  are  on  opposite  sides  of  the  stock  price. 
Although  standard  and  partial  barrier  options  use  flat  barriers,  double  barrier  options 
where  the  barriers  may  be  curved  are  treated  in  [8].  In  order  to  price  such  an  option,  the 
density  must  reflect  that  the  stock  has  not  hit  either  of  the  barriers.  They  use  an  extension 
of  Levy’s  density  for  Brownian  motion  confined  to  an  interval  [a,b]  (i.e.  a Brownian 
motion  for  which  a <m,  < M,  <b).  This  density  is  based  on  repeated  applications  of  the 
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reflection  principle  of  Brownian  motion  and  is  written  as  an  infinite  sum  [10],  The  density 
is 


{x+2j{b-a)f 


It 


-e 


(x-v2j{b-a)-2b) 
2t 


2 N 


dx 


(2.5) 


where  A - < m,  < M,  < b,B,  g 1 c [a,  i]} , Their  results  indicate  that  the  convergence 

of  this  series  is  sufficiently  rapid  to  make  it  amenable  to  numerical  implementations. 
Despite  the  generality  of  their  results  with  curved  boundaries,  we  will  confine  ourselves  to 
double  barrier  options  in  which  the  barriers  are  flat.  In  what  follows,  we  will  refer  to  these 
types  of  barrier  options  as  tunnel  options. 

Just  like  other  barrier  options,  double  barrier  options  may  be  used  to  speculate  on 
the  perceived  volatility  for  an  asset.  There  are  many  different  ways  to  compute  volatility, 
none  of  which  are  universally  accepted  [5].  Consider  an  investor  who  expects  the  volatility 
to  be  less  from  implied  volatility  of  a given  call  option.  He  could  sell  the  call  option  and 
purchase  the  same  option  with  an  appropriate  double  barrier.  If  he  is  right,  then  both 
options  will  be  worth  the  same  at  expiry  and  will  thus  leave  him  with  no  obligations.  He 
will  profit  by  the  amount  equal  to  the  difference  in  the  options’  prices.  Double  barrier 
options  have  been  developed  and  traded  by  several  Tokyo  banks  including  Fuji  Bank, 
Sanwa  Bank  and  Nippon  Credit  Bank  [8]. 


CHAPTERS 

SOME  NEW  RESULTS  ON  BROWNIAN  MOTION 

This  chapter  presents  some  new  results  on  several  functionals  of  Brownian  motion. 
These  results  will  be  of  great  use  in  chapter  4.  The  corollaries  following  the  main  theorems 
provide  densities  for  a Brownian  motion  at  time  T confined  to  a tunnel  during  a fixed 
period.  The  fixed  period  will  either  start  at  time  zero  and  end  at  some  time  t before  time 
T or  will  begin  at  some  time  t and  end  at  time  T . We  use  several  lemmas  to  conceal 
tedious  computations  and  obstacles  in  the  proofs  of  main  theorems. 

Lemma  3.1  Fix  the  values  x,,X2  andy,  and  define  the  function  h{c,d)  by 

\2  A 


h{c,d)  = 


1 


/ f f exp 


cy- 


(x  + df  (y  - x) 


2t 


2{T-t) 


dydx 


(3.1) 


Then 


h{c,d)  = exp 


( c^T 

cd 

j 

P 

\ 

P-yx 
It  ’ Vr 


\ V { ^-^2  P-yx'^ 

^ /t  ’ 4f 


(3.2) 


where  (•,  •)  denotes  the  bivariate  normal  distribution  function  with  correlation 


coefficient  p - , a = ct-d  and  f -cT  -d . 


(x-^-df  —Y 

Proof:  Let  f{x,y)  = cy ^ . We  first  factor  out  the  quantity 


2/  2(7-0 


2(7-0 


from / so  that  we  have 
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f(x,y)  = 


-T 


2(7-0 


2(7-0  , 2(T-t){y-xf  ^ 2{T-t){x  + df 

T ^ T 2(7-0  7 2t 


• (3.3) 


Expansion,  cancellation,  grouping  like  terms  and  simplification  yields 


-7 


2(7  - 1) 


x^+2d 


^-7 

Tj 


2xy  ^ y -2c(T-t)y 
7 T 


2t 


-T7-  (3-4) 


We  would  now  like  to  express  the  quantity 


r 


X + 2d 


1--  X 


V 


2x;^  ^ y -2c(J-t)y 
T T 


(3.5) 


in  the  form 


{x-af  2(x-a)(y-P)  , {y-pf  , 

h h Y . 

t T T 


(3.6) 


In  order  to  solve  for  the  constants  a,  P and  y we  first  equate  the  coefficients  of  x and  y 
terms  from  (3.5)  and  (3.6).  This  leads  to  the  system  of  equations 


2d 


u 

t 


2a  2P 
t T 


-2c{T-t)^  2p  ^2a 
7 ~ T T 


(3.7) 


Solving  this  system  yields 


a -ct-d 
p^cT-d' 


(3.8) 


We  next  set  the  constants  from  both  expressions  equal  and  obtain 


a"  2ap  p^ 
-y  = — + 7— . 

/ 7 7 


(3.9) 
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So  we  may  now  rewrite  (3.4)  as 


-T 


2{T~t) 


{x-af  2{x-a\y-P)  {y-pf  , _ 

h h V 

t T T 


(3.10) 
2t  ’ 


where  a, >5 and/  are  defined  as  in  (3.8)  and  (3.9).  We  take  the  / term  outside  of  the 
brackets  so  that 


2{T  - 1) 


(,x-a)  2(x-a)(y-p)  ^ jy-pf 

t T T 


-T  (f 

H r 

2{T-ty  2t 


(3.11) 


-T  dP 

We  now  simplify  the  term  — — — -/ . Substituting  the  values  for  a and  p from  (3.8) 

2(  i /)  2/ 

into  (3.9)  for  y , the  above  expression  becomes 


-T  P 

-y-—  = 


J f 0^/3  /?2  A „2 


2(7-0  2/  2(7-0 

7 


ay  _ 2^  ^ p^ 


c 

Tt 


2(7-0 


{ct-df  2{ct-d){cT-d)  jcT-d) 
t T 7 


2 A 


2t 


(3.12) 


After  expanding  and  simplifying  (3.12)  we  have 


-7 


d^  PT 


-y = cd . 

2(7-0  2 


(3.13) 


Substituting  (3.13)  into  (3.11)  yields 


fix,y)  = 


-7 


2(7-0 


{x-af  2{x-a){y-p)  ^ {y~pf 
t T T 


Pt  , 

^ cd 


Set  a --  \ft,  a - \[f  and /?=./—  . Then  we  have  — = 

V7 


(3.14) 


— and  P - ^ 


2(7-0  T 


so  that  (3.14)  becomes 
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f{x,y)  = 


-1 


(x-g)  2p{x  - a){y  - p)  ^ {y~pf 
O'."  o-,or^  o-/ 


c'r  . 

H cd . 


(3.15) 


Recall  that 


+00 

f 

ix+df 

J J 

•ti  y, 

^y 

2t 

2(T-t)  J 

dydx . Thus  we 


have 


1 .J  +00 


6Y 


6c 


27ryjt(T-t) 


nyh. 


JJ 


{x-af  _2p{x-aXy-p)_^  jy-pf 

a/  CT^y  CT„^ 


(3.16) 


dydx 


R y 

Making  the  substitutions  x = , y = — — 


o. 


and  switching  the  limits  of  integration 


yields  (3.2).  □ 

Lemma  3.2  The  density 


k{x)  = 


1 


\llKt 


;=-cc 


exp 

(x  + 2j{b  - a)f 

-exp 

{x  + 2j{b-a)-2bfy 

2t 

2t 

. 

is  uniformly 


convergent  on  the  interval  [a,  b] . 
Proof:  If  a series  of  the  form 


Z exp 


oo 


(x  + 2j{b  - a)f 
2t 


(3.17) 


converges  uniformly  on  the  interval  [a,b],  then  a series  of  the  form 

+00 

Z e^P 


j — ® 


{x  + 2j{b-a)-2by 
2t 


(3.18) 
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will  also  converge  uniformly  on  the  same  interval.  To  see  this,  first  recall  that  a <0<b  . 
Then  observe  that 


{x  + 2j(b-a)f  ^ (x  + 2U  + l)(b-a)-  2bf 
2t 


> -- 


2t 


for  large  positive  values  of  j and 


{x  + 2j{b-a)f  ^ {x  + 2j{b-a)-2bf 
2t  ^ 2t 

for  large  negative  values  of  j.  So  the  convergence  of  (3.17)  ensures  that  (3.18)  will  also 
converge.  Thus  we  need  only  show  that  a series  of  the  form  (3.17)  converges  uniformly  on 
the  interval  [a,Z>] . However,  this  is  equivalent  to  the  one  sided  series 


exp 


(x  + 2j{_b  - a)f 
2t 


(3.19) 


converging  uniformly  on  the  same  interval.  Now  we  choose  an  integer  N large  enough 
that  j >N  implies  x + 2 j(b  -a)>  0.  We  then  observe  that  for  j > N we  have 


exp 


(x  + 2jjb  - a)f 
2t 


exp 


{a^2j{b-a)f 

2t 


Vxe[a,Z>].  (3.20) 


Now  the  series 


exp 


(g  + 2j{b  - a)f 
2t 


(3.21) 


clearly  converges.  If  we  define 


r 


M„  = sup 

jce(a,fc] 


exp 

7=1 


(x  + 2y(Z>-g)) 
2t 


(3.22) 


then  we  have 
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f 


lim  = lim  sup 

„-*+C0  "-*■«>  ;ce[a,i] 


2;  exp 


]-n 


(o  + 2y(6  - a)) 


2 lA 


2/ 


Jy 


< lim 

n—^+oo 


^exp 


y-n 


{a  + 2j{b  - a)) 
It 


2 lA 


(3.23) 


= 0 


The  inequality  in  (3.23)  is  due  to  (3.20).  The  limit  goes  to  zero  since  the  series  converges. 
(3.23)  is  a necessary  and  sufficient  condition  for  the  uniform  convergence  of  (3. 19)  [11]. 
This  proves  the  lemma.  □ 

Theorem  3.3  Let  A = {a<m,  <M,  <b,Bj  > >'j } where  t <T . Then  we  have 


E 


e^Bj  . j 


A 


X (Kx,cj)-h(x,d,)) 

y=-co 


where  h{-,  ■)  is  given  by  (3.2),  Cj  = 2 j(b  - a)  and  dj  = 2 j{b  -a) -2b. 


(3.24) 


Proof:  Let  ={a  <m,  < M , < b]  and  4 = [Bj  > } so  that  A = A^r^A2  - Then  we 

have 


e^^  1 


■ M 


= E 


A,  1 


(3.25) 


Using  a simple  property  of  expectations,  we  take  an  inner  expectation  conditioned  on 
so  that  (3.25)  becomes 


E 


E 


(3.26) 


We  now  note  that  the  random  variable  1^^  is  .7^ -measurable.  Thus  we  can  take  1^ 
outside  of  the  inner  conditional  expectation  so  that  (3.26)  becomes 


Ia  E 


,XBt 


(3.27) 
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Next  we  rewrite  (3.27)  using  the  shift  operator  6,  and  the  Markov  property  of  Brownian 
motion.  This  yields 


= E 


1.  E 


e^-  \^  \-e, 


e^^-‘  ■ 1 


(3.28) 


where  A2  = {Bj_,  > y'j } and  E^  [•]  denotes  the  expectation  for  the  process  started  at  z 
instead  of  the  origin.  We  may  compute  the  inner  expectation  as 


which  is  a function  of  B, . So  we  define  the  Sanction  / by 

fix)  = E 

and  write  (3.28)  as 


2{T-t) 


dy,  (3.29) 


j +00 


2(7-0 


dy  (3.30) 


^A,m) 


Recall  that  /4,  = [a  <m,  < M,  < Z>} . So  we  compute  (3.3 1)  as 


where 


k{x)  = 


flnt 


(x  + 2j{b  - a)f 

[x  + 2j{b-a)-2b)'^ 

exp 

2t 

-exp 

2t 

(3.31) 


(3.32) 


(3.33) 


Set  Cj  = 2 j{b  - a)  and  dj  - 2 j{b  -a)-2b . Then  we  may  express  (3.32)  as 
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J 


1 +00 


exp 


It 


exp 


2t 


f{x)dx.  (3.34) 


Substituting  /(x)  back  into  (3.34)  and  pulling  the  constants  to  the  front  of  the  summation 
we  obtain 


1 


1 


2n^t{T-t)  jTL 


exp 


2t 


- exp 


(x  + c/j' 


2t 


\j  exp(^>')  • exp 


f ( \2  ^ 


(3.35) 


dydx 


2{T-t) 

where  / = [>^i,oo) . Since  /(x)  is  bounded  on  the  interval  [a,b],  we  use  lemma  3.2  to 
interchange  the  infinite  limit  with  the  outer  integral.  We  may  also  combine  the  integrands 


to  obtain 


J CO 


=11. 


27i^t(T  -t) 


exp 


2t 


-exp 


2 


2t 


r r 

exp(A>')  • exp 


V 


ilzA. 

HT-t) 


2 


dydx 


JJ 


(3.36) 


Distributing  the  exponentials  on  the  right  yields 


I 


b 

\ 

!L 

a 

exp 

2t  2(7-0 


exp 


2t 


2(7-0 

We  now  have  a doubly  infinite  series  where  the  terms  are  of  the  form 


dydx 


(3.37) 


h{X,Cj)-h(A,dj) 


(3.38) 
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where  /?(•,  ) is  defined  as  in  (3.2).  We  use  lemma  3.1  to  complete  the  proof.  □ 

+00 

Corollary  3.4  P{a  <m,  <M,  <b,Bj  > y^)=  (h{Q,Cj)-h{Q,dj))  for  t < T . 


Proof:  This  is  immediate  from  theorem  3.3.  □ 

Again  we  fix  the  constants  Xj, Xj  and>',  and  define  a new  function  h(c,d)  by 


1 -*2  •*'1 
Hc,d)  = , f j 


exp 


c>^- 


2 > 


2t  2(7 -0 


dydx.  (3.39) 


Remark;  This  function  differs  from  h{c,d)  only  by  the  limits  of  integration  on  they 
variable. 

Lemma  3.5  The  function  h{c,d)  may  be  expressed  as 
h(c,d)  = exp 


( 7\ 

y 

cd 

1 2 J 

P 

\ 

where  p 
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= a=c,~ 


7?  ' 77 

d and  f -cT  -d  . 


^2-«  f^i-«  y\-P^ 

■'’p 


y 


4t  ’ Vr 


y 


(3.40) 


Proof:  This  proof  is  the  same  as  the  proof  of  lemma  3.1  except  for  making  the  last 
substitutions 


_ x-a  . y-f 

X = and  y - - — — 

<y.  O',, 


(3.41) 


instead  of 


- a-x  . p-y 

X = and  y = — — 

a.. 


(3.42) 


□ 


Lemma  3.6  Define  the  function  g(c,d)  by 
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2 -tj  >2 

= 1 f f exp 

2n4I(T^){{ 


2 A 


cy- 


{x-^df  (>'-y) 
It  2{T-t) 


dydx.  (3.43) 


Then 

g{c,d)  = ex^ 


( c'^T 

r 

f f 

cd 

J 

LI 

[ % 

x^-a  y^-p 

VF  ’ VF  , 


-N„ 


N„ 


X2-a  y,-p 

, VF  ’ Vf 


^ r 
■N„ 


Xj-a  y2~P 

>/F  ’ VF  , 

^1  -g  Ti -P 

’ VF  , 


(3.44) 


where  (•,  •)  denotes  the  bivariate  normal  distribution  function  with  correlation 
coefficient  p - , a = ct-d  and  f = cT  -d . 

Proof:  This  proof  is  immediate  since  g{c,d)  is  the  difference  of  two  functions  having  the 

, ■‘2  >’2 

form  h{c,d)  or  h{c,d)  In  particular,  any  integral  of  the  form  1 1 {^(fydx  may  be 


*1  yi 


^2  >2 


written  as  J"  J ( )dydx  -j  j"  (^-^dydx  so  that  the  right  hand  side  is  the  difference  of  two 


functions  of  the  form  h(c,d) . □ 
Theorem  3.7  Define  w{c,d)  by 


w{c,d)  = exp 


TV, 


cd 


\ 


''xj  +d-a  f Xj  +d-a  y^-p^ 


V 


VF  ’ VF 


y 


V 


VF  ’ VF 


y 


X2  +d-a  2izZl_Ar  ( +d-a  y,~P 


VF  ’ VF 


VF  ’ VF 


(3.45) 
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andsetB-{a<B^<b'^ue[t,T],Bj&l]where  t<T  and  I = [y^,y^~](z[a,b].  Then 
we  have 


r 1 

E[e^ -\sy 


(3.46) 


j^—oo 


where  Cj  = 2 7'(Z)  - a)  and  dj  - 2 y(Z>  -a) -2b. 

Proof:  Just  as  in  theorem  3.3  we  use  basic  properties  of  conditional  expectation  and  the 
Markov  property  of  Brownian  motion  to  obtain 


= £ 
= E 


(3.47) 


where  B - {a  < B^  < 6 Vw  g [0, 7 - t],Bj_,  e /} . We  may  compute  the  inner  expectation 


as 


(3.48) 


where 


k{x) 


1 +00 

(x  + 2j(b  - a)f 

- exp 

{x  + 2jib-a)-2bfl 

^27T(T-t)j±L 

2(7-0 

2(7-0 

(3.49) 


Define  the  function  / (x)  by 

f{x)  = E''  e^^-'  = ^jk{y  - x)-  Qx^[Xy)dy  . 


(3.50) 


We  may  then  compute  the  entire  expectation  from  (3.47)  as 
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£[£''[£-‘«'-.1j]]  = £[/(5,)] 

We  now  substitute  the  function  f{x)  into  (3.51)  to  obtain 

1 * 

1 1/  ~ (^y)  ^ 

and  the  density  A:(x)  into  (3.52)  which  becomes 


v’2^ 


dx 


(3.51) 


v'20 


dx 


(3.52) 


II; 


yj2n{T  - 1) 
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exp 


2 A 


2(7’-/) 
&x\)[Xy)dy  ■ exp 


exp 


2(7’-/) 


2 ^ 


(3.53) 


v’2^ 


dx 


Now  / (x)  is  bounded  so  we  may  use  the  uniform  convergence  of  the  series  to  move  the 
summation  outside  of  the  integral  signs.  After  simplifying  (3.53)  becomes 

[y-x  + c^f  x^ 


,._oo  27iyjt{T 
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f 

Jl, 

a 

exp 

V 

Ay-- 


2(7’-/)  2/ 


exp 


Ay- 


2 

(y-x  + dj) 
2(7’-/)  2/ 


(3.54) 


dydx 


Define 


w{c,d)  = 


2nyjt{T  -t) 


JC,  Vj 

I Jexp 

■^i  y> 


x^  (y-x  + d) 
cy ^ 


2 ^ 


2/  2(7’  - /) 


dydx 


(3.55) 


and  recall  that 


g(c,^7)  = 
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yi 

(x  + ^)^ 

2n^tiT- 

1)  J J ®"^P 

■'i  >’i 

2/ 

2(7’-/)  J 

^7va!)c.  (3.56) 
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Notice  that  we  could  translate  the  x variable  in  w(c,d)  and  rewrite  (3.55)  as 


2 ^2-‘i  yi 

w(c,d)  = ■■  f fexp 


cy- 


{x  + df  (>^-x) 


2 A 


2t  2(7-0 


cfydx . (3.57) 


We  appeal  to  lemma  3.4  to  see  that  w(c,d)  may  be  written  as 


w{c,d)  = exp 


- cd 


-d-a  -p 

/t  ’ Vr 
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y 


^1  -d-a  y^ -p 

yft  ’ Vr  jj 


V 


X2-d-a  yi~P]_j^  (x^-d-a  y^-p^ 

■‘''p 


' Vr 


V 


.Tt  ’ ^/7 


y 


(3.58) 


where  p - Jy  , a = ct-d  and  P -cT -d  . But  then  w(c,d)  = w(c,d)  so  (3.54) 


becomes 


+<X) 


(3.59) 


This  completes  the  proof.  □ 


Corollary  3.^  P [a  < <b  \/u  e I)  - ^ (vr(0,  ) - ^(0,  d^ ))  for  t < 7 

y=-QO 

and  I = [yx,y2]^W,b]. 

Proof:  This  is  immediate  from  theorem  3.6.  □ 


CHAPTER  4 

PARTIAL  TUNNEL  OPTIONS 
4,1  Description 


A partial  tunnel  option  is  precisely  what  its  name  indicates.  An  option  of  this  type 
does  not  monitor  the  barrier  for  the  entire  life  of  the  option.  Just  as  with  a partial  barrier 
option,  the  monitoring  period  either  starts  at  time  zero  and  ends  at  a certain  time  before 
expiry  or  starts  at  a certain  non-initial  time  /j  and  ends  at  expiry.  We  distinguish  these 
two  cases  as  Type  I and  Type  II  options,  respectively.  This  chapter  will  give  explicit 
pricing  formulas  for  call  and  put  options  of  both  Type  I and  Type  II  in  terms  of 
cumulative  bivariate  distribution  functions. 

4.2  Removing  the  Drift 


Recall  from  chapter  1 that  the  price  of  an  option  is  simply  the  expectation  of  its 
discounted  payoff  with  respect  to  the  martingale  measure  P . The  payoffs  from  each  of 
the  partial  tunnel  options  depend  on  the  distribution  of  the  stock  price  process 


S,=Se 


aW,A 


r-- 


(4.1) 


where  W,  is  a Brownian  motion  under  P . All  of  the  densities  we  have  derived  in  chapter 
3 and  most  that  are  found  in  the  literature  are  based  on  Brownian  motion  with  no  drift 
term.  However,  the  Brownian  motion  that  appears  in  the  stock  price  process  is  always 
coupled  with  a drift  term.  Observe  that 
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C7fV,+ 


In  particular,  the  process  we  are  dealing  with  is  Brownian  motion  with  drift  — 

a 


(4.2) 


r - 


In  order  to  use  densities  for  different  functionals  of  Brownian  motion  we  must  remove 
the  drift  term.  This  can  readily  be  done  using  the  Cameron-Martin-Girsanov  theorem  as 
we  did  in  chapter  1.  We  implicitly  define  the  probability  measure  Q , under  which  the 
process 


W^=W,+- 

a 


r- 


t 


(4.3) 


is  a Brownian  motion,  via  the  Radon-Nikodym  derivative 


dP  ^ 


^ ^ 2 


= exp|^//(lT7- 
= exp  nWj  — fi^T 


(4.4) 


where 


1 

G 


r - 


(4.5) 


Observe  that  the  quantity  n defined  in  (4.5)  is  not  the  original  drift  for  the  stock  price 
process  as  in  chapter  1 . Here  fx  is  the  drift  of  the  stock  price  process  under  the 
martingale  measure  P . Now  the  arbitrage  price  for  an  option  is  found  by  taking  the 
expectation  under  the  martingale  measure  P . The  expectation  under  the  measure  P is 
not  equal  to  the  expectation  under  the  measure  Q in  general.  However,  the  Radon- 
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Nikodym  derivative  serves  to  relate  the  two  expectations.  In  particular,  let  X . Then 
we  have 


[X]  = E9 


(4.6) 


Thus  the  arbitrage  price  for  an  option  may  be  computed  as  an  expectation  under  Q of  the 
discounted  payoff  multiplied  by  the  Radon-Nikodym  derivative  in  (4.4).  In  this  manner 
we  will  be  able  to  use  the  densities  for  Brownian  motion  with  no  drift. 

4.3  Partial  Tunnel  Options  of  Type  I 
4 3.1  Partial  Tunnel  Call  Option  of  Type  I 

A partial  tunnel  call  option  of  type  I (PTCO-I)  has  the  payoff  function 


where  k is  the  strike  price  and  the  barriers  are  a and  b . Thus  the  arbitrage  price  of  a 
PTCO-I  is  given  by 


= E^ 


dP 


(4.8) 


where  the  second  equality  is  from  (4.6).  Under  the  new  measure  Q we  have 


S,  = S e 
= Se 


aWi+ 


(4.9) 


aiV, 


Define  A-{a<  S^<b  Vm  e [0,  /,  ] } . Using  (4.9)  we  rewrite  A as 


A=[a<  S^<b  Vmg[0,/,]} 
= [a<  W^<b  Vm€[0,/,]} 


(4.10) 
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In 


where  a = — 


In 


s 


and  b - — . Substituting  from  (4.4)  and  A from  (4.10)  we 
(T  dQ 


may  express  the  price  from  (4.8)  as 


exp 


f ~ \ ^ 


^-rT 


^a<W^ib  VMe[0,/|]|  ^ 


Lemma  4.1  The  price  of  the  PTCO-I  is  given  by 


J,P+(T)Wj 


where  k = 


f k\ 

In  - 

-All 


-k-£f 


V«e[0.r,],tFj.6[;t,oo)j 


(4.11) 


(4.12) 


Proof:  Since 


{.Sr -k>0}={Wj.  >yt} 

= \Wr^[k,^)] 


exp 


pW,-^p^T 
2 y 
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Vue[0,/,]| 

Next  we  break  this  up  into  two  separate  expectations  so  that 


(4.13) 


• (4.14) 
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pW, 
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pW, 

\ 
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\a<W^<b  Vtte[0,/,].lf're[/fc,<»)|‘ 


2 


L - - - - J 


(4.15) 


Replacing  Sj  with  the  process  in  (4.9)  and  factoring  out  constants  yields 
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S-E^ 


•1 


Vwe[0,/,].Pr7.e[i^,oo)| 


-k-EP 


This  proves  the  lemma.  □ 

Theorem  4.2  The  price  of  the  PTCO-I  is  given  by 


(4.16) 


+00  —rT——iTT 

2 \^s\h{p  + a,Cj)-h(,p  + a,dj))-k-[h{p,Cj)-h(,p,dj))'^  (4.17) 

/— 00 

where 


h{c,d)  = exp 


f c^T  ,] 

N 

a-b  /3-k 

-N 
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•Q 

1 

1 2 J 

p 

•'''p 

(4.18) 


Cj  - 2{b  -d),  dj  = 2{b  -d)-2b , p = Jy  , a = ctj  -d  and  f-cT-d 


Proof:  From  lemma  4. 1 we  have 


V 


•1 


-k-EP 


a<Wy^b  Vu6[0,/|],^rj.e[^,«)| 

VMe[0,/,],#j.e[^,oo)| 


(4.19) 


Applying  theorem  3.3  to  each  of  the  expectations  completes  the  proof  □ 

The  price  of  the  PTCO-I,  as  well  as  the  price  of  each  of  the  other  partial  tunnel 
options,  is  written  out  explicitly  in  the  appendix. 

4.3.2  Partial  Tunnel  Put  Option  of  Type  I 

Computations  similar  to  those  above  can  be  used  to  show  that  the  arbitrage  price 
Vf  for  the  partial  tunnel  put  option  of  type  I (PTPO-I)  can  be  expressed  as 
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-rT-\n^T  ^ 


Vf=e 


k-E^ 


^a<fV^^b  Vi/e[0,/,],fFj.e(-oo,;fc)| 


-S-E^ 


,(M+(tWt 


• 1 


a^fP’u^b  Vi/e[0,r,],#j.e(-oo^)| 


(4.20) 


jy 


Theorem  4.4  The  price  Vf  of  the  PTPO-I  is  given  by 


+«>  —rT——jET  r / -. 

^ \j^\Kf^^^j)-h{p,dj)^-S-(h{p  + <j,Cj)-h{p  + a,dj)^.  (4.21) 

where 


h{c,d)  = exp 

Cj^2{b-a),d^=2{b 

Proof:  We  have 


cd 


^ b-a  k - 


-N^ 


( ~ r ^ ^ 
a-a  k-  p 


(4.22) 


-a)-2b , p = , a = ct^-d  and  p = cT-d 


k-E^ 


V 


-S-E^ 


jaSfF„<6  Vi/e[0./,],fFj.e(^Jt]} 

^(p+a)fVj.  j 


Using  the  same  methods  as  in  proof  of  theorem  3.3  we  can  see  that 

+00 


(4.23) 


(4.24) 


where  A - ^ < m,^  - ^r,  <b,Bj-e  (-°o,>'i]|.  Using  (4.24)  to  evaluate  (4.23)  completes 


the  proof  □ 
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Remark:  The  only  differences  between  the  valuation  of  the  PTCO-I  and  the  PTPO-I  are 
the  payoffs  and  the  limits  of  integration.  This  is  what  motivated  the  introduction  of  the 
function  h{c,d). 


4.4  Partial  Tunnel  Options:  Type  II 
4 4. 1 Partial  Tunnel  Call  Option  of  Type  II 

The  same  type  of  computations  can  be  used  to  show  that  the  arbitrage  price 
for  the  partial  tunnel  call  option  of  type  II  (PTCO-II)  can  be  expressed  as 


S-E^ 


V«6[r, j'],irj.e[fc,cxD)| 


-k-E^ 


eP^T  . I 


(4.25) 
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Theorem  4.5  The  price  of  the  PTCO-II  is  given  by 


^ -rT-^fpT 


= Z ^ ^ \s  [Hp  + ^.Cj)~w{p  + G,dj))-k-{w{p,Cj)-w{p,dj))'^{A26) 


where 


w{c,d)  = exp 
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b +d-a  b- p 
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a + d-a  b- p 

E ’~JT, 


b+d-a  S~P 

ylh  ’~VF 


.X  .(4.27) 


-N„ 


d+d-a  5- p 
yjr,'y[f 


5 = max(a,^),  c^  -2(b-d),  dj  = 2(b -d)-2b  , a = ct^  -d  and  p = cT-d  . 
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Proof:  We  have 


-k-EP 


•1 


Vm6[/, j’],fTj.e[it,ao)| 


(4.28) 


We  do  not  know  whether  k <a  or  a <k  . This  is  the  reason  S = max(a,k)  is 
introduced.  We  may  then  rewrite  (4.28)  with  ^ yielding 

jM+cr)fVj. 


\ 


^a^fV^<b  Vwe[/,j'],tFj.e[ty,oo)| 


-k£P 


(4.29) 


We  then  apply  theorem  3.6  to  evaluate  the  expectations  which  yields  (4.26).  This 
completes  the  proof.  □ 

4 4.2  Partial  Tunnel  Put  Option  of  Type  II 


The  same  type  of  computations  can  be  used  to  show  that  the  arbitrage  price 
for  the  partial  tunnel  put  option  of  type  II  (PTPO-II)  can  be  expressed  as 


p 

II 


-rT--u^T  ( 

2^  k-E‘^ 


•1 


-S-E^ 


^ai.Wu'^b  VMe[r,,7'],#j.e(a^]| 


(4.30) 


Theorem  4.6  The  price  of  the  PTPO-II  is  given  by 


where  w{c,d)  is  as  defined  as 
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w{c,d)  - exp 


-cd 


b +d -a  y/  - P 
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VL  V 


<K  ’ ^ 


-N. 


^ a + d -a  y/  - p^ 

•F, 


^ b +d -a  a- p^ 

4h 


-N. 


f 

a+d -a  a- p 


,(4.32) 


y/  = min(b,k),  Cj  =2(b-a),  d^  =2(b-a)-2b,  a = ct^ -d  and  p = cT-d . 
Proof:  We  have 


V 


|a<#„SA  VMe[/|,r],fFj.6(a,^]| 


•1 


VMe[/,  j'],fF^e(a,/t]| 


(4.33) 


jy 


Again  we  do  not  know  whether  k <b  or  b <k . So 'wt  introduce  y/  = min(i,  k)  and 
rewrite  (4.33)  as 


V 


k-E^ 


.1 


Vi/e[r, 


-S-E^^ 


(4.34) 


We  then  apply  theorem  3.6  to  evaluate  the  expectations  which  yields  (4.3 1).  This 
completes  the  proof.  □ 


CHAPTER  5 

ANALYSIS  OF  THE  PARTIAL  TUNNEL  OPTIONS  PRICING  FORMULAS 

The  first  part  of  this  chapter  shows  how  the  partial  tunnel  options  pricing  formulas 
generalize  those  of  the  existing  options  pricing  formulas.  Two  methods  are  used,  both  of 
which  are  straightforward  limiting  procedures.  In  particular,  we  take  the  limit  of  a partial 
tunnel  pricing  formula  as  the  upper  and  lower  barriers  go  to  infinity  and  zero,  respectively. 
We  also  fix  the  upper  and  lower  barriers  and  take  the  limit  as  the  monitoring  time  goes  to 
zero.  In  the  case  of  Type  I partial  tunnel  options,  both  of  these  limits  turn  out  to  be  the 
pricing  formulas  for  the  standard  options.  However,  a subtlety  arises  in  the  latter  limit  for 
partial  tunnel  options  of  Type  II.  We  first  revert  back  to  the  original  expressions  for  the 
options  prices,  as  expectations  of  payoffs.  This  will  greatly  simplify  the  task  as  we  take 
limits.  Second,  we  take  the  limit  of  the  new  pricing  formula  given  in  chapter  4 as  an 
infinite  series.  Although  this  requires  much  more  detail  and  cumbersome  analysis,  we  will 
see  the  Black- Scholes  formula  emerge  as  the  limit  of  the  PTCO-I  price  as  the  monitoring 
period  /j  approaches  zero.  The  second  part  of  the  chapter  provides  numerical  results  and 
explains  the  patterns  found  in  the  results. 
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5, 1 Partial  Tunnel  Options  Pricing  Formulas  as  Extensions  of  Existing  Formulas 
5.11  Limits  of  Partial  Tunnel  Options  Prices  as  Expectations 

Recall  that  the  arbitrage  price  of  a PTCO-I  is  given  by 


^-rT 


(5.1) 


Clearly  we  have 


< (5^  - Arf 


Moreover, 


and 


(5.4) 


Since  these  limits  are  increasing  and 


^{a<S,<b  V«e[0,/,]}  (‘^7’  ^ 0 


(5.5) 


we  may  use  the  Monotone  Convergence  Theorem  (twice)  to  obtain 


lim 

a-*0 

b-*+OD 


v«e[0,r,]}(‘^7'  {St  k) 


(5.6) 


That  is,  the  limit  of  the  price  of  a PTCO-I  option  as  the  barriers  go  off  to  oo  and  0 is  the 
price  of  a standard  call  option. 

Remark:  The  first  application  of  the  Monotone  Convergence  Theorem  shows  that  the 
partial  tunnel  option  price  converges  to  that  of  a partial  barrier  down  and  out  call  as  the 
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upper  barrier  goes  to  cc . Since  the  limits  can  be  taken  in  either  order,  we  also  see  that 
the  partial  tunnel  option  price  converges  to  that  of  a partial  barrier  up  and  out  call  as 
the  lower  barrier  goes  to  0. 

Another  way  to  see  (5.6)  is  to  take  the  limit  as  the  length  of  the  monitoring  period,  /j , 
goes  to  zero.  In  particular,  we  have 

=('^7-"^)  • (5-7) 

Since  this  is  also  an  increasing  limit  and  we  have  (5.5),  we  may  again  apply  the  Monotone 
Convergence  Theorem  to  obtain  (5.6). 

Remark:  Clearly  we  could  repeat  the  same  analysis  for  the  PTPO-I.  However,  in  the 
case  of  the  Type  11  partial  tunnel  options,  taking  the  limit  as  t^  approaches  T does  not 
converge  to  the  price  of  a standard  option.  In  this  case  the  payoffs  are  always  bounded. 
This  is  due  to  the  fact  that  all  of  the  terms  in  the  limit  are  uniformly  bounded.  No  matter 
how  close  /,  is  to  T , the  barriers  still  exist  at  time  T and  therefore  impose  a bound  on 
the  payoffs  for  call  and  put  options.  They  do  however  converge  to  the  prices  for  call-like 
and  put-like  options  whose  payoffs  are  given  by  X = 1^^^^  [Sj  - k)*  and 

5 . 1 .2  Limits  of  Partial  Tunnel  Options  Formulas 

In  this  section  we  will  focus  on  the  PTCO-I.  In  particular,  we  will  take  the  limit  of 
the  pricing  formula  obtained  in  chapter  4 as  the  length  of  the  monitoring  period  /j 
approaches  zero.  The  formula  is  expressed  as  a doubly  infinite  series.  We  will  see  that  all 
of  the  terms  except  for  the  central  term,  the  term  corresponding  to  j = 0,  will  approach 
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zero  and  therefore  drop  out  from  the  formula.  The  resulting  limit  is  the  classic  Black- 
Scholes  formula  for  a call  option.  Note  that  we  use  the  uniform  convergence  from  chapter 
3 to  interchange  the  infinite  sum  with  the  limit. 

The  formula  for  the  PTCO-I  found  in  the  appendix  is  given  by 
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In  — 

where  - 2(b  -d),  dj  = 2{b  -d)-2b  and  p = .1—  . Recall  that  d = — — and 

''  r a 


In- 

- C . 

b Since  a < S and  b>  S we  have  d <0  and  b >0  . First  we  consider  the 

o 


following  quantities 
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-dj-b 

for  7 0 . In  particular,  we  will  show  that  each  is  either  positive  or  negative  according  to 

whether  j is  positive  or  negative.  We  substitute  = 2{b  - a)  and  dj  = 2(b  -a) -2b 
into  (5.9)  so  that  we  have 


-Cj - a = -2j{b-d)-d, 

-Cj  -b  = -2j{b-a)-b, 

-dj  - a = -2  j(b  -a)  + 2b  - a 
-dj  - b = -2  j{b  -d)  + 2b  -b. 


(5.10) 


Simplifying  yields 


-Cj  - d = {2j-\)a-2jb, 
-Cj-b  = 2jd  - {2j  + ])b, 
-dj-a  = (2j-\)a-(2j-2)b 


(5.11) 


-dj-b  = 2ja-i2j-2)b. 

Observe  that  if  j >0  then  all  of  the  quantities  in  (5. 1 1)  are  negative.  On  the  other  hand,  if 
j <0  then  all  of  the  quantities  are  positive.  These  observations  should  be  evident  from  the 


fact  that  a <0  and  b > 0 . Now  consider  the  difference 


^ + (j)T  - Cj 


^ -b  {f^+a)T -Cj -ic^ 


(5.12) 
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For  j > 0,  the  first  argument  in  each  bivariate  normal  distribution  fimction  is  negative  for 
sufficiently  small  values  of  . As  /j  tends  to  zero  each  of  these  arguments  goes  to  -oo , 
Thus  each  distribution  function,  and  hence  the  entire  expression,  will  go  to  zero.  For 
j <0,  the  first  arguments  are  both  positive  and  therefore  tend  to  -foo . Each  bivariate 
distribution  function  then  converges  to  a (univariate)  normal  distribution.  Since  the  second 
arguments  in  each  bivariate  distribution  function  are  the  same,  the  difference  approaches 
zero.  The  same  argument  shows  that  each  of  the  three  other  difference  expressions  go  to 
zero  as  well.  This  leaves  us  with  the  term  corresponding  to  j = 0.  This  term  simplifies  to 
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Eight  bivariate  normal  distribution  functions  appear  in  this  term.  We  use  the  same  analysis 
as  we  did  above.  Recall  that  a < 0 and  b>0  . Note  also  that  2b-a>0 . As  tj  goes  to 
zero,  the  first  and  fifth  bivariate  distribution  functions  converge  to  a (univariate)  normal. 
The  second  and  sixth  each  go  to  zero.  The  four  remaining  distribution  functions  each  go 
to  1 and  therefore  cancel  each  other  out.  Thus  the  first  term  is  simply 
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Replacing  /^  = — 
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r - 
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we  may  simplify  (5.14)  to  recover  the  classic 


Black-Scholes  formula. 


5.2  Numerical  Results 


In  this  section  we  examine  some  features  of  the  PTCO-I  and  PTCO-II.  Naturally, 
the  PTPO-I  and  PTPO-II  will  share  many  of  these  features,  so  they  will  be  left  out.  For 
each  of  the  partial  tunnel  call  options,  we  will  numerically  illustrate  the  analytic  results  of 
section  5 1 We  will  also  examine  the  effects  that  the  volatility  and  the  length  of  the 
monitoring  period  have  on  the  price  of  each  option. 

5.2.1  The  Partial  Tunnel  Call  Option:  Type  I 

Consider  a PTCO-I  where  we  have  the  following  parameters;  initial  stock  price 
6’  = 55,  volatility  cr  = .2,  expiry  T = 1,  monitoring  time  /,  = .5 , strike  price  k -6S, 
lower  barrier  a = 40,  upper  barrier  ^ = 80  and  interest  rate  r = .06 . We  will  first  show 
that  as  the  upper  and  lower  barriers  of  this  option  approach  infinity  and  zero,  respectively. 
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the  price  of  the  option  will  approach  that  of  the  standard  call  option  with  initial  stock  price 
6’  = 55 , volatility  a = .2 , expiry  T = \,  strike  price  k = 65  interest  rate  r = .06 . The 
appropriate  price  for  this  standard  call  option  is  $2.166 . We  price  the  same  PTCO-I  with 
upper  barriers  ranging  from  80  to  120  and  lower  barriers  ranging  from  40  to  5.  The  results 
are  presented  in  figure  5.1.  Observe  that  the  price  is  the  same  as  the  Black-Scholes  price 
when  the  lower  and  upper  barriers  have  been  moved  to  40  and  120 , respectively. 


We  now  consider  the  case  above  but  let  the  length  of  the  monitoring  time  go  to 

zero  In  order  to  illustrate  the  effect  that  this  has  on  the  option  price,  we  will  use  the  same 
parameters  as  above  except  that  the  upper  and  lower  barriers  will  be  70  and  50 , 
respectively.  We  observe  the  prices  with  /,  ranging  from  0 to  1 . The  results  are  displayed 
in  figure  5.2.  Observe  once  again  that  the  price  approaches  the  standard  Black-Scholes 
price  as  approaches  zero. 


(40,  120,  2.166) 


pnce 


Figure  5.1:  Varying  upper  and  lower  barriers  for  a PTCO-I 
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Figure  5.2:  Varying  the  length  of  the  monitoring  period  for  a PTCO-I 


We  now  observe  the  affect  that  the  volatility  a and  the  length  of  the  monitoring 
period  t,  have  on  the  price.  In  the  case  of  the  standard  call  option,  increasing  the  volatility 
increases  the  price  of  the  option.  However,  in  the  case  of  the  partial  tunnel  options, 
increasing  the  volatility  only  increases  the  price  up  to  a point.  After  this  point,  the 
probability  of  hitting  a barrier  and  having  a zero  payoff  is  so  high  that  increasing  the 
volatility  only  decreases  the  chance  of  not  hitting  a barrier  and  having  a positive  payoff. 
Thus  increasing  the  volatility  after  this  point  only  serves  to  decrease  the  price.  This 
phenomenon  is  present  in  both  Type  I and  Type  II  partial  tunnel  options. 

The  length  of  the  monitoring  period  is  also  a factor  in  the  price.  For  partial  tunnel 
options  of  Type  I the  monitoring  period  begins  at  time  zero  and  ends  at  some  time  . 

Thus  the  length  of  the  monitoring  period  is  simply  /, . Since  increasing  only  increases 
the  chance  of  getting  a zero  payoff,  we  see  that  the  price  will  decrease  as  we  increase  . 
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Here  we  take  the  case  where  the  intial  stock  price  is  = 20 , cr  = . 1 8 , r = .06,  7 = 2 , the 
upper  barrier  is  45,  the  lower  barrier  is  15  and  the  strike  price  isyt  = 25 . Figure  5.3  plots 
the  values  of  the  prices  as  the  volatility  ranges  from  .05  to  .95  and  the  monitoring  time 
ranges  from  0 to  2. 


5 
4 

price  3 
2 
1 

Figure  5.3:  Varying  cr  and  /j  foraPTCO-I 

Notice  that  as  the  length  of  the  monitoring  period  gets  closer  to  zero,  the  price 
appears  to  increase  as  <j  increases.  This  occurs  as  a result  of  the  PTCO-I  converging  to  a 
regular  call  option  where  the  price  always  increases  as  volatility  increases.  This 
phenomenon  will  not  occur  in  options  of  Type  II. 

5. 2. 2 The  Partial  Tunnel  Call  Option:  Type  II 

Consider  a PTCO-II  where  we  have  the  following  parameters:  initial  stock  price 
‘S’  = 55 , volatility  a - .2,  expiry  T -\,  monitoring  time  /,  = .5 , strike  price  A:  = 65 , 
lower  barrier  a = 40 , upper  barrier  = 80  and  interest  rate  r = .06 . We  show  that  as  the 
upper  and  lower  barriers  of  this  option  approach  infinity  and  zero,  respectively,  the  price 


monitoring  period 
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of  the  option  will  approach  that  of  the  standard  call  option  with  the  same  parameters.  We 
will  not  be  taking  the  limit  as  the  length  of  the  monitoring  time  goes  to  zero  for  the 
reasons  mentioned  in  the  above  remark.  Recall  that  the  appropriate  price  for  this  standard 
call  option  is  $2. 166 . We  price  the  same  PTCO-II  with  upper  barriers  ranging  from  80  to 
125  and  lower  barriers  ranging  from  55  to  25.  The  results  are  presented  in  figure  5.4. 
Observe  that  the  price  is  the  same  (to  at  least  three  digits  of  accuracy)  as  the  Black- 
Scholes  price  when  the  lower  and  upper  barriers  have  been  moved  to  25  and  125, 
respectively. 


We  again  observe  the  effect  that  the  volatility  cr  and  the  length  of  the  monitoring 
period  have  on  the  price  options.  Now  for  options  of  Type  II,  the  monitoring  period 
begins  at  /,  and  ends  at  time  T . Therefore  the  length  of  the  monitoring  period  in  this  case 
is  r - /j . As  mentioned  above,  increasing  the  volatility  only  increases  the  price  up  to  a 
point  and  then  the  price  begins  to  decrease.  We  now  look  at  the  the  case  where  the  intial 
stock  price  is  ^ = 55 , a = ,2 , r = .06,  T = 2 , the  upper  barrier  is  100,  the  lower  barrier  is 


(25,125,2.166) 


1.5: 
price  ^ : 


Figure  5.4:  Varying  upper  and  lower  barriers  for  a PTCO-II 
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30  and  the  strike  price  is  ^ = 65 . The  volatility  ranges  from  .05  to  .30  and  the  length  of 
the  monitoring  period  ranges  from  0 to  2.  These  results  are  presented  in  figure  5.5  below. 
Note  that  the  axis  for  the  monitoring  period  represents  the  length  of  the  monitoring  period 
r - 1,  which  was  /j  in  the  case  of  the  Type  I option  above. 


Figure  5.5:  Varying  a and  T-ti  for  aPTCO-II 


Observe  the  plot  where  the  length  of  the  monitoring  period  T gets  closer  to 
zero.  As  a increases,  the  price  still  only  increases  to  a point  and  then  decreases.  This 
phenomenon  is  different  from  that  of  the  PTCO-I  where  the  price  appeared  to  strictly 
increase  as  a increased  for  small  values  of  . 


CHAPTER  6 
THE  GET-OUT  OPTION 

61  Description 

The  get-out  option  is  based  on  a strategy  that  will  get  one  out  of  the  market  when 
a certain  trigger  is  initiated.  Consider  a bullish  investor  who,  for  some  reason,  cannot 
monitor  the  market  and  wishes  to  implement  a strategy  that  would  exit  him  from  further 
exposure  if  one  of  his  holdings  were  to  go  below  a certain  level.  The  get-out  option 
depends  upon  two  underlying  securities  and  its  payoff  is  determined  by  whether  or  not 
these  securities  go  below  a certain  level  during  the  lifetime  of  the  option.  There  are  three 
possible  occurrences; 

1 . Stock  1 hits  a lower  barrier  a and  triggers  the  payoff  of  an  exercised  call  on 
Stock  2 before  time  T. 

2.  Stock  2 hits  a lower  barrier  b and  triggers  the  payoff  of  an  exercised  call  on 
Stock  1 before  time  T. 

3.  Neither  stock  hits  its  respective  barrier  before  the  expiry  time  r and  the 
option  returns  the  payoffs  of  call  options  on  both  stocks  at  time  T. 

This  option  is  different  from  the  options  of  previous  chapters  in  that  it  depends  on 
two  underlying  stocks.  As  a result,  the  Black-Scholes  setup  from  chapter  1 will  no  longer 
suffice.  We  introduce  a generalization  of  this  model  that  will  enable  us  to  price  this  new 
option  appropriately. 

6.2  The  Multidimensional  Black-Scholes  Setup 

6 21  The  Model 

We  must  now  use  the  multidimensional  (2  dimensional)  Black-Scholes  Model  [9]. 
In  this  model  we  have  the  filtered  probability  space  (Q,  T,  P)  in  which  the  processes 
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dS,=S,{u,dt  + i:,-dW,) 
dB,  = VfBjdt 


(6.1) 


are  defined.  Here  S,  is  a (2x1)  vector  of  stock  price  processes,  u,  is  a (2x1)  vector  of 
drifts,  I,  is  the  volatility  (2x2)  matrix  and  W,  is  a multidimensional  (2x1)  Brownian 
motion  under  the  probability  measure  P . We  also  have  the  natural  filtration  that  is 
the  filtration  generated  by  the  2 dimensional  process  S, . 


Definition  6.1  A process  y for  which  r,  • 1 - m,  = Z,  • y,  is  called  the  market  price  for  risk. 


The  existence  of  this  process  in  conjunction  with  Girsanov’s  Theorem  gives  rise 
to  a martingale  measure  for  our  model.  Recall  that  the  martingale  measure  is  simply  the 
measure  under  which  the  discounted  stock  price  process  is  a martingale.  Since  the  stock 
price  process  is  2 dimensional,  this  implies  that  each  component  of  the  process  e~’''S,  is  a 
martingale  For  simplicity  we  consider  the  case  where  u,  = u,  Z,  = Z is  nonsingular  and 
r,  = /• . In  particular,  the  drift,  volatility  and  interest  rate  do  not  depend  on  time.  Under 
these  assumptions  it  is  clear  that  the  process  y indeed  exists  and  is  given  by 
y - IT^r -\-u) . This  model  together  with  these  assumptions  is  referred  to  as  the 
classic  Black-Scholes  model. 

6.2.2  The  Martingale  Measure 

As  before,  we  begin  with  the  martingale  measure  P . Under  this  measure  we  have 


'si'" 

(6.2) 
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for  i=l,2  where  Z'  denotes  the  i***  row  of  the  matrix  Z . In  particular,  the  discounted  stock 
price  process  is  a martingale  under  P . The  solution  to  this  stochastic  differential  equation 
is  given  by 


f 

f 2 2 > 

\ 

= S'  exp 

Z'  • W,  + 

_ 0-,l  +CT/2 

t 

1 ^ ) 

y 

(6.3) 


when  Z is  given  by 


V*^21 


(6.4) 


We  may  write  out  (6.3)  as 

S',  = S'  exp 


f 

f 

2 2 > 

aX  + 

0-,i  +0-,-2 

t 

2 J 

V 

K 

y 

(6.5) 


6.2  3 The  Payoff 

The  discounted  payoff  X for  the  get-out  option  can  be  expressed  as  the  sum  of 
three  exclusive  payoffs  that  correspond  to  the  three  possible  outcomes  set  forth  in  5. 1 . Set 

r,  = inf  ^5  > 0 : S]  = ^|,  ^2  = inf  > 0 : S^  = and  define  the  following  sets 

accordingly: 

A^  ={r,  <T2;  r,  <7’} 

^2={h<h'y  (6.6) 

4 ={r,  >T,  T2>T} 

Clearly  we  have  P [A^  A2  <J  Aj)  = \ and  o = 0 for  j . So  the  discounted 


payoff  X may  be  expressed  as 
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(6.7) 


Our  goal  now  is  to  price  this  option.  Just  as  before,  we  need  only  take  the  expectation  of 
X under  the  martingale  measure  P . 

6.3  Pricing 

6.3, 1 The  Independent  Case 

This  is  the  simpler  of  the  two  cases.  For  this  case  we  assume  that  the  two  stock 
price  processes  are  independent.  Consequently,  the  non-diagonal  entries  in  the  volatility 
matrix  Z must  be  zero.  Otherwise  the  processes  would  both  contain  the  same  nonzero 
Brownian  motion  terms  and  therefore  be  correlated.  If  we  denote  the  diagonal  entries  of 
the  volatility  matrix  by  ct,  = cr^^  and  CTj  = CTjj  then  we  may  write  the  stock  price 
processes  as 


S;  exp(o-,JF/ 


(6.8) 


where  //,  = r - 


for  i-1,2.  Recall  from  (6.6)  and  (6.7)  that  we  will  be  taking  the 


expectation  of 


(6.9) 


where 
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A ={^1  <^2',  ^\<T] 
A ={^2  <^I>  ^2  <^j 

4 = {r,  >7;  T2  >r} 


(6.10) 


We  begin  with  the  first  term.  We  wish  to  compute 

E 


“ 

f ^ 

/ \ + \ 

'4 

e ' 

1 

1 

(6.11) 


under  the  martingale  measure  P . Recall  that  the  process  Sl  is  defined  as 

Si  = expfo-jlf;^  + //2«) 

Define  the  processes  by 

+ n-u  for  i=l,2. 


(6.12) 


Then  we  have 


SI  =S‘  exp(A'')  for  i=  1,2. 


(6.13) 


(6.14) 


Now  write  the  set  ^ as  4 = A'r^A"  where  4'  = {r,  <72}  and  4'=  {rj  < Tj.  If  we 
have  the  density  P(tj  g du,  < T2,XI^  e (fy)  then  we  will  be  able  to  express  (6. 1 1)  as  a 

Lebesgue  integral.  First  note  that  4'  = > In  ^ | where  mf  is  the  running 

minimum  for  the  process  Xf . Thus  we  may  write 

P(t,  g du,  r,  < T2  , X^  &dy)^  P(r,  g du,  m\  >h,X]  Gdy)  (6.1 5) 


where  ^ = In 


VO  y 


Using  methods  in  [12]  and  the  independence  we  have 


P(r,  G du,  > b,  X]^  ^dy)=  P(r,  g du,  ml  >b,XlG  dy) 

= P(r,  G du)  ■ P(ml  >b,XlG  dy) 


(6.16) 


53 


The  densities  on  the  right  hand  side  are  found  in  [7],  The  density  p(^Xl  &dx,ml>  is 
given  by 


P(xUdx,ml>b)  = 

1 


TUJ^U 


exp 


2<7^U 


exp 


V y 


exp 


Icr^u 


dx 


■ (6.17) 


Recall  that  r,  - inf  ^5  > 0 : S]  = urj  We  substitute  5"]  = exp^X]  ^ and  rewrite  this  as 


T^  = inf  ^5  > 0 : 5"]  = a| 

= inf  > 0 : exp^Zj)  = a| 

= inf  ^5  > 0 : Xl  = a| 

= inf|5>0:  IT/  +^s  = a[ 


(6.18) 


In 


where  a = 


v«y 


Then  we  have  the  density 


P(r,  e du)  = 


\a-^u] 

\d\ 

1 J 

1 ■ V 

2u 

du 


(6.19) 


So  we  may  express  (6. 1 1)  as 


T +00 

= j J e"™  ) P(r,  G du,  r,  < Tj, e dy)  (6.20) 

0 —00 
T +00 

= I j ) P(r,  G du)P(^X^  edy,ml> 

0 —00 
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We  remove  the  (•)'^  notation  by  changing  the  lower  limit  of  integration  of  y to 


k-,  = In 


so  that  the  quantity  inside  of  (■)'"  above  is  strictly  positive.  Then  (6.20) 


becomes 


T +00 


J j e ™ (S^e^  - ^2)^  P(Ti  € du)P{^Xl  e dy,ml  > b)  . (6.21) 

0 

Notice  that  the  second  expectation  is  the  same  as  the  first  with  the  roles  of  stock  1 and 
stock  2 reversed.  Thus  we  can  write  it  as 

T +00 

1 J P(t2  g du)P(^Xl  ecfy,ml>  a)  (6.22) 


0 t, 


where  A:,  = In 


fir  ^ 


and  <5  = In  I 

S' J 


Lastly,  we  need  to  compute  the  value  of  the  third  expectation 


= E 


|(*^r  ^1 ) + {Sf  A^2  ) 


(6.23) 


We  now  distribute  the  indicator  function  to  each  of  the  terms  in  the  parentheses  to 
compute  each  expectation.  The  first  term  is 


■ E 


% 


{s'r-k,y 


(6.24) 


Recall  that  /I3  = {r,  > T;  r2  > T} . We  may  write  ^3  as  n A”  where 


^3  - {^1  > and  A^  - {^"2  > T).  We  break  up  the  indicator  function  accordingly  and 


express  (6.24)  as 
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-rT 


■E 


^ 


(■^r  ~^i) 


(6.25) 


We  may  factor  this  into  two  separate  expectations  using  the  independence  so  that  (6.25) 
becomes 


Note  that  we  have 


e 


-rT 


(6.26) 


E 


= p{Ai) 

= P(z,>T) 


(6.27) 


■Jl 


T 


exp 


v2  ^ 


b-^S 


2s 


ds 


\ y 

The  last  equality  is  given  in  [7].  Thus  (6.26)  is  simply  the  price  of  a down  and  out  barrier 
option  multiplied  by  the  probability  in  (6.27).  This  formula  is  given  in  [9].  The  second 
term  in  the  expectation  (6.23)  is  the  same  with  the  roles  reversed. 

Remark:  Notice  that  the  above  analysis  did  not  use  the  change  of  measure  employed  in 
earlier  chapters  to  remove  the  drift.  The  reason  for  this  was  that  the  densities  used 
already  accounted for  the  processes  with  drift. 


6 3 2 The  Dependent  Case 

Just  as  before,  under  the  martingale  measure/*  the  stock  price  processes  are 


S[  - S'  exp 


r - 


2 , 2 ^ ^ 

q~,i  +^,2 


for  i=L2 


(6.28) 


where  JT/  and  fV,^  are  standard  Brownian  motions  . Moreover,  we  are  still  trying  to 


compute  the  expectation  (under  P)  of  the  discounted  payoff  N where 
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^ = 1 


A 


-n. 


K ^2) 


and 


+\ 

■(42 

-*.)1 

(6.29) 

+^A, 
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•)) 

A 

= {2’l 

<T2; 

z,  <Tj 

A 

“ {^2 

<2-1; 

r,<r). 

(6.30) 

A 

^={h 

>T- 

IV 

^a} 

and  T 

2 = inf  > 0 : 

. This  time  we  will 

remove  the  drift  with  a change  of  measure.  Notice  that  we  must  remove  the  drift  term 
from  each  of  the  processes  in  (6.28).  We  will  change  to  a measure  Q so  that 


where  W,^  and  are  standard  Brownian  motions  under  Q . However,  the  drift 
constants  ju^  and  jU2  must  be  chosen  so  that  they  satisfy 

- Mi‘)  + o-,2  - M2*) 

= o-,,#;  + a,2W,^  - o-,,//,/  - 
= ^,1^/  + - a, 2^2)^ 


(6.31) 


(6.32) 


r - 


^2  , _2  ^ 


for  i-1,2.  In  particular,  and  ^2  satisfy  the  equations 


^12/^2  ~ ^ 


2 2 

<^11  <^12 


2 2 
CT21  +<^22 


(6.33) 


+ ^22/^2  ~ ^ 


2 
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Solving  these  equations  for  //,  and  //j  yields 


2r  {022  cT\2 ) ^22  ^12 ) ^12  i^2i  ■*“  *^22 ) 

i i i L 

2(cTi,Cr22  ~0'j2(T2j) 


(6.34) 


M2  = 


_ 2/-  (gn  - 0-21 ) - 1 (q-2^  + 0-22  ) + 0-21  (o-U  + 0-12  ) 


2(<Jjjg22  <^12*^21 ) 


So  we  define  a new  measure  Q via  the  Radon-Nikodym  derivative 


dP 

— = exp 
dQ 


V 


(6.35) 


J 


Under  the  measure  Q the  stock  price  processes  are 


S',  = S'  exp(c7„fr/  + ct,2^/^)  for  i^l.2 


(6.36) 


where  W,^  and  W,^  are  standard  Brownian  motions  under  Q . We  may  now  compute  the 
expectation  of  X from  (6.29)  as 


[X]  = EP 


^X 

dQ 


(6.37) 


Recall  that  r,  = inf  ^5  > 0 : S]  = « j We  may  rewrite  this  as 

r,  = inf  > 0 : 5’’  exp {^<y^^W^  + a, 2^/ ) “ 
= inf  > 0 : CFnWj  + ct,2^/  = 


(6.38) 


where  a = In 


. Similarly,  we  may  write  Tj  = inf  |5  > 0 : g2i^i  +^22^^  = where 


b = \n 


r h\ 
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We  will  make  two  substitutions,  a translation  and  a rotation,  so  that  we  may  use 
the  joint  densities  found  in  [6]  to  compute  the  expectation  in  (6.37).  Note  first  that  the 
stopping  times  rj  and  correspond  geometrically  to  the  first  time  a two  dimensional 


Brownian  motion  W,  - 


, starting  at  the  origin,  first  hits  either  one  of  the  lines 

Lj  : o'iiX  + o'i2y 

L^.  CT2iX  + a22y  = b 


(6.39) 


Note  that  a <0  and  b <0  since  a <S^  and  b<  . So  the  intersection  of  these  lines  lies 


in  the  third  quadrant.  The  point  of  intersection  p = 


Uo, 


of  these  two  lines  is  given  by 


^22^  (Tj2^ 

^11^22  ~ *^12*^21 


(6.40) 


(T t tb  o' 'y\Cl 

^0  = ^ 

CTj,(T22  CTi2<^21 


Denote  the  angle  between  the  two  lines  as  a (note  that  this  is  the  obtuse  angle  containing 
the  origin  since  both  lines  have  negative  slopes). 


Figure  6. 1 : Original  two  lines  in  the  plane 
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The  densities  in  [6]  are  for  a process  that  starts  at  specified  point  in  the  plane  and  then 
hits  one  of  two  lines.  One  of  these  lines  is  the  x-axis.  The  other  line  goes  through  the 
origin  and  has  a negative  slope.  The  point  at  which  the  process  starts  must  be  located  in 
the  wedge  that  includes  the  first  quadrant.  Thus  we  will  need  to  relocate  the  intersection 
of  our  two  lines  to  the  origin  and  then  rotate  the  plane  so  that  one  of  the  lines  coincides 
with  the  x-axis.  In  this  manner  we  will  be  able  to  use  the  densities  from  [6]. 

We  now  make  the  first  substitution 

^.=W,-p  (6.41) 

so  that  Xj  is  a Brownian  motion  starting  at  -p  . This  translation  moves  the  intersection 

of  the  two  lines  to  the  origin.  The  slopes  of  line  1 and  line  2 are  /w,  = 

<Ti2 

and»J2  = respectively.  We  may  suppose  without  loss  of  generality  that  w,  > W2 . 

022 

Note  that  we  could  simply  switch  the  roles  of  stock  1 and  stock  2 to  accomplish  this. 
However,  the  assumption  that  the  slopes  are  not  equal  may  not  be  as  easy  to  see.  If  the 
slopes  are  equal  then  the  stocks  are  perfectly  correlated.  Since  this  case  is  of  no  interest 
we  may  make  the  assumption  that  the  slopes  are  indeed  not  equal. 


Figure  6.2:  Translated  plane 
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The  next  substitution  rotates  the  plane  counter  clockwise  so  that  line  1 coincides 


with  the  x-axis.  Since  the  slope  of  line  1 is  /Wi  = 11- , the  angle  that  it  makes  between 


'12 


itself  and  the  positive  x-axis  is 


6 = tan 


-1 


’-'11 

V ^12  J 


(6.42) 


Note  that  we  have  ^ < 0 since  — ll-  < 0 . Note  also  that  the  angle  between  line  2 and  the 


'12 


x-axis  is  also  negative.  Thus  the  angle  is  given  by 


a - n - tan' 


CT21 

-i-tan  * 

f-— 1 

V *^22  J 

1 *^12  y 

(6.43) 


n 


We  have  a>-^  since  both  lines  have  negative  slopes.  Define  a new  process  Y,  by 


y.=R-o-x, 


(6.44) 


where  R_n  is  defined  as 


R-e- 


^ cos  9 sin^'' 
-sin^  cos^ 


(6.45) 


The  matrix  R_(^  will  rotate  the  plane  counterclockwise  by  an  angle  of  -9  > 0 . Thus  7,  is 
a two  dimensional  Brownian  motion  starting  at  the  point  R_g  ■(-p)  = -R_gp . Recall  that 
the  original  process  W,  starts  at  the  origin  and  the  intersection  of  the  two  lines  is  in  the 
third  quadrant.  Thus  the  translated  and  the  rotated  process  Y,  starts  in  either  the  first  or 
second  quadrant.  Moreover,  line  1 becomes  the  x-axis  and  line  2 creates  an  obtuse  wedge 
that  includes  the  first  quadrant  and  a portion  of  the  second  quadrant.  This  comes  from  the 
fact  that  the  original  lines  both  have  non-positive  slopes.  Note  that  one  of  the  slopes,  but 


61 


not  both,  may  be  zero.  In  this  case  the  angle  between  line  1 and  the  positive  x-axis  is 
zero.  Thus  the  rotation  matrix  above  turns  out  to  be  the  identity  matrix  in  this  case. 


Figure  6.3:  Translated  and  rotated  plane 


Now  the  process  7,  may  be  expressed  as 

Y,=R.e-X, 

= R-9-[W,-p)  . (6.46) 

= R-e-W, -R_ep 


Solving  backward  for  W,  yields 


W,=RZl{Y,+R,p) 

= R-ey,^p 


(6.47) 


We  now  look  at  what  the  individual  processes  and  look  like  in  terms  of 


the  process  7, . So  we  write  out  (6  47)  as 


W,=RZl-Y,^p 


-sin^'' 

fv'  ^ 

•*  J 

+ 

fx  ^ 

cos^  ^ 

7^ 

Vi  ^ 

.yoy 

^cos0-  7/  - sin  ^ • 7,^  + ^ 

^sin  ^ • 7/  + cos  9-Y,^ 


(6.48) 
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We  define  the  functions  / and  g by 

f(x,y)  = xcos^- j/-sin^  + Xo 


and 


^(x,  >^)  = X • sin  ^ + >^  • cos^  + >^o . 


(6.49) 


(6.50) 


Then  we  have 


w,= 

(6.52) 


(6.51) 

This  will  help  us  write  the  increasingly  large  expressions  more  compactly.  Now  the 
stopping  time  r,  defined  as 

= inf  ^5  > 0 : ctjjIF/  + 

may  be  expressed  as 

Tj  = inf  |5  > 0 : CTjj  (cos^  • Y'  - sin  ^ • 7/  + Xq) 

+<Ji2  (sin  0 • 7/  + cos  ^ • 7/  + j = a| 

A simple  calculation  shows  that 

0-11^0 +0-1270  =«  (6.54) 

Thus  we  may  write  (6.53)  as 

r,  = inf{5>0:  a,,  (cos6»  ■ 7/  - sin  (9  • 7/ ) + <7,2  (sin  6>-7/  +cos(9-7/)  = oj 


(6.53) 


= inf|5>0.  (cr,,  cos  +<7,2  sin  6>)  7/  + (cr, 2 cos  (9 -a,,  sin  (9)  7/  = o| 


• (6.55) 


Moreover,  d = tan' 


V 0-12 


SO  we  may  compute  the  trigonometric  functions  sin  0 and 


cos^  as 
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sin^  = 


cosff  = 


-cr. 


II 


+0-12 


'12 


yj^n  +0-] 


(6.56) 


Substituting  these  into  (6.55)  yields 

T,  = inf  > 0 : • 7/  = o|  (6.57) 

Clearly  the  quantity  + a^2  is  strictly  positive,  thus 

r,  = inf  {5  > 0 : 7/  = o} . (6.58) 

This  is  precisely  what  we  should  expect.  Recall  that  r,  represents  the  first  time  that 
process  7,  hits  line  1 . However,  line  1 has  been  the  translated  and  rotated  so  that  it  is  the 
x-axis.  Thus  r,  represents  the  first  time  that  the  process  7,  hits  the  x-axis  (i.e.  the  first 
time  that  7/  = 0 ). 


Our  goal  is  to  compute  the  price 
this  out  explicitly  as 


dQ 


X, 


of  the  get-out  option.  We  first  write 


dQ 


X 


dP 


f 


dQ 


i 


•(■??, -*2)* 


(6.59) 


+ ] 


/ / 1 \ / 

\ + 

((sj-*,)  +(j 

T ~ ^7) 

)))] 

The  first  term  in  this  expectation  is 


dQ 


-n. 


'K  ^2) 


(6.60) 
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Under  the  measure  Q we  have 

S\  = 5' exp^CTjjfr/  +<t,2^/)  for  i=  1,2 
and 

dP 


(6.61) 


dQ 


= exp 


(6.62) 


where  W,  and  are  standard  Brownian  motions.  We  now  make  the  substitutions  for 
IF/  and  W,^  from  (6.51).  Thus  (6.60)  becomes 


exp 


-n. 


(6.63) 


But  r,  = inf  > 0 : = o|  so  = 0 . Replacing  this  above  yields 


1^"  exp  (o-2,  ■ / , o)  + 0-22  • ^ (y/j , o)j  - ^2 ) 


'A 


~f~y\ 
e ' 


(6.64) 


We  use  some  simple  properties  of  expectations  to  take  an  inner  expectation  conditioned 
on  the  sigma  algebra  . This  yields 


E^ 


E 

f 

exp 

V 

.2  , ..2  ^ 


^.f[Yl,Y,^)  + M2g(yrJT)-^^^^T 

exp  (cT2i  ■ / (y/j  , o)  + CT22  • g (y/j  , 0 jj  - ^2 1 


Er  r 

./  aTi 


(6.65) 
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We  take  out  the  measurable  portion  inside  of  the  inner  expectation.  All  that  is  left 

is  the  Radon-Nikodym  derivative,  which  is  a martingale.  Recall  that 
A = {^1  < ^2.  ^1  < 7’} . In  particular,  <T  on  this  set  so  T a Tj  = r,  on  this  set.  Using 
the  Optional  Sampling  theorem  [7]  for  the  process  7,  and  the  bounded  stopping  time 
r A Tj , (6.65)  becomes 


■ exp  ,oj  + a22-g  , o))  - j 


exp 


^ J 


= E^ 


exp 


■ exp  (<7„  ■ / , o)  + <7„  ■ g , o))  - J 


f^/{Yr,.Yf^)^M2g{Yl^.YT])- 


2 2 \ 


(6.66) 


We  again  observe  that  7/^  = 0 and  substitute  this  into  (6.66)  to  obtain 


E^ 


exp 


exp  (ct2,  ■ / , o)  + aj2  • g (7^^ , o)j  - j 


(6.67) 


Observe  that  the  random  variable  inside  of  the  expectation  depends  only  on  r,  and  7^^ . 
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So  we  define  a function  H (■,•)  by 


H{u,a)  = e (5^exp((T2i  ■ f {a,0)  + 0^2  ■ g{a,Q))  - 

^ 2 , 2 ^ 

exp  f{a,0)^  ^g{a,Qi)-^  9 

I ^ 


(6.68) 


Note  that  the  (-)^  notation  has  been  removed.  We  may  do  this  as  long  as  we  make  sure 

that  the  quantity  inside  ( )^  is  greater  than  or  equal  to  zero.  The  following  equivalent 
statements  illustrate  how  to  ensure  this. 

1)  5^exp(CT2,  /(a[,0)  + cr22-g(a,0))-^2  ^0 


2)  0-2,  /(a,0)  + o-22-^(a,0)>A2  =ln 


(h.] 

U'J 


3)  CT21  (acos^  + Xo)  + CT22(asin^  + ^o)>^2  = 


4)  a2i  ■ a cos  9 + CT22  • a sin  ^ + (ct2,Xq  + CT22>'o  ) - ^2 


5)  (T21  * Cl  cos  ^ "1“  ^'2,'^  * ^ Sin  0 b ^ h-2 


6)  a(cT2i  cos^ + CT22  sin^)  > A:2  - i 


7)  a 


f ^ 

^12^21  ~ *^11^22 


> k.,  - b 


o\  ^ I 2 ^ 

8)  a> yjcTu  + C7i2 

(Ti2<T2i  O'n‘^22 

Statements  1-4  and  6 come  from  substituting  and  rearranging,  5 comes  from  a calculation 
similar  to  that  of  (6.54)  and  the  coefficient  on  or  in  6 is  positive  so  we  may  divide  both 
sides  by  it  to  get  to  7. 
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Define  the  constant 


Tj  = max 


k^-b 


yCTijCTji  CTjjCr22 


+0-12,0 


(6.69) 


If  ^2  > ^ , then  we  will  only  need  to  compute  the  expectation  of  H j over  values 
of  which  are  greater  than  ^2  If  ^2  - ^ > then  we  will  compute  the  expectation  of 

H j over  all  values  of  Yj^ , i.e.  values  greater  than  or  equal  to  zero.  So  the 

constant  is  a lower  bound  for  the  values  of  Y^^  for  which  we  must  take  the 
expectation.  Thus  we  may  write  (6.67)  as 


h i.  . 


But  the  joint  density  for  the  pair  ) when  r,  < T2  given  in  [6]  is 

&da,T^  <r2j 


(6.70) 


n 


2 , 2 ^ 

a +/-n  ' 


— exp 

a au  \ 2u 


Z 

;=0L 


7 -Sin 


V cc  J 


V ‘ 7 


(6.71) 


where  p - 


''tq  cos6q^ 

^ To  sin  ^0  7 


, a - 7T  -tan 


-1 


^21 

+ tan~‘ 

^11 

V ^22  7 

1*^12  7 

and  Ip  is  the  modified 


Bessel  function  of  order  P . Recall  that 


+00 


/ \2i 

' X ' 


k\T{p  + k + \) 

p+00  _ , 

where  T(a)  = e f dt . Recall  that  a is  simply  the  angle  between  the  lines  in 
which  the  process  Y,  starts.  Thus  we  may  express  (6.70)  as 


(6.72) 
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T +00 


\ \ H {u,a)p(x^  e duJl^  eda,r^<  Tj) 
r +00  r 


2 , 2 ^ 
a +/-n  ' 


0 /: 


a au 


(6.73) 


z 

y-0 


7 -sin 


. a / 


^a/- ^ 


V ' y 


dadu 


This  representation  allows  for  its  precise  value  to  be  calculated  via  standard  numerical 
integration  techniques  as  opposed  to  simulation. 

This  is  only  the  first  term  in  the  expectation.  Fortunately,  the  same 
transformations  enable  us  to  use  the  densities  in  [6]  to  write  out  the  expectations  as 
Lebesgue  integrals.  The  second  expectation  we  must  compute  is 


(6.74) 


We  again  use  the  functions  / and  g defined  in  (6.49)  and  (6.50)  so  that  we  have 


W = 


fwy 

^y.'.y.^). 

(6.75) 


We  may  then  write  (6.74)  out  as 


exp 


[s'  exp(a„ 


(6.76) 


We  use  the  same  method  as  in  (6.66)  so  that  (6.76)  becomes 


exp 


2,2  ^ 
Ml  +P2  , 

~T~  \ 
^ y 


f ^ 1 / 

/ \ / \\ 

\+V 

e 2 .|^iexp|cr„ 

(6.77) 
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Recall  that  = inf  ^5  > 0 : + ctjjPF/  = i j . So 

Tj  = inf  ^5  > 0 : (72iW^  + 0’22^/  = 

= inf  {5  > 0 : • / (n* , 7/ ) + • g {y^  J^)  = b] 

= inf  > 0 : CT21  (7/  cos^-7/  sin  ^ +0-22  (j",*  sin^  +7/  cos^  + Xq)  = 

= inf  ^5  > 0 : (T21  (7/  cos  ^-7/  sin  ^ + Xq  ) + <722  (7/  sin^  +7/  cos^  +Xq)  = 

= inf  |5  > 0 : (<T2i  cos0  + a 22  sin  6)Yl  (6.78) 

+ (o-22  cos^  - CT21  sin  e)Y^  + (0-21X0  + a22J^o)  = b] 

= inf  |5  > 0 : (a-21  cos 6 + O22  sin  &) 7/  + (a22  cos 6 - (J21  sin 9) Yf  +b  = 

= inf  > 0 ; (a2i  cos^  + a 22  sin  6)Yl  + [a 22  cos 9-  0-21  sin  9)Y^  = o| 

= inf|5>0:  7;  =^11^21+^120-22  .^4 
I O'!  1 0-22  “CTj  20-21  J 

The  coefficient  of  the  7/  term  is  precisely  the  reciprocal  of  the  slope  of  line  2 after  the 
transformation.  Define  another  constant 


^ _ 0^110^21  +0^120-22 

0"llO-22  ~0-i2CT2i 


(6.79) 


so  that 

r2  = inf |5>0:  7/  = . yA 

I O11CT22  - cTi20-2i  J (6.80) 

= inf{5>0:  7;  =/w-7/} 


Then  we  have 


(6.81) 


Now  we  may  rewrite  (6.77)  as 
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exp 


exp  (a„  • / (m  • , i"r2  ) + 0-12  • ^ (^  • >^4  ’ ))  " ^ ) 


-,  (6.82) 


so  that  the  random  variable  we  are  taking  the  expectation  of  depends  only  on  Tj  and  7^^ 
Define  a function  M (• , ■)  by 


2 , ..2  ^ 


M(w,a)  = exp  //,  ■ f{m-a,a)  + fi^-g{m  a,a)-^  u 


[e  exp(o-„  •/("»  a,a)  + ai2  .g(/w  a,a))-^i)) 


(6.83) 


Then  we  must  again  find  the  lower  bound  rj  for  the  random  variable  so  that  we  may 
omit  the  (-)^  notation.  Computations  similar  to  those  for  /j  yield 


r 


X2  = max 


k^-a 


I 2 2” 

+cr,2 


,0 


(6.84) 


So  we  may  express  (6.82)  as 


(6.85) 


Thus  we  only  need  the  joint  distribution  for  the  pair  ^X2,  Y^^  j given  r2  < r, . But  this  is 


given  in  [6]  by 
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p(t^  sdu,Y^^  Gda,T2  <r,j 
;r-cos^a-^j 


^22 
a +Tq  -cos 


2 exp 

a au 


2u  ■ cos^ 


1 

a-‘ 

yi) 

/2)  y/ 

(6.86) 


I 

2-0 


-1^  • y-  sin 


^ jn_^ 
V « y 


ar. 


COS 


where  9^  - a - 6^.  Note:  If  q is  the  point  we  get  by  reflecting  p = 


cos  6^ 

V^bsin^o  j 


about  the 


f - 

line  _y  = tan  — ■ x , then  q - 

V 2 J 


/•q  cos  9q 

V'bsin^ 


. The  cos^a  - factor  appears  as  a result  of 


projecting  the  process  Y,  onto  the  y-axis  to  obtain  the  random  variable  Y^^  . 

Remark:  The  formula  in  [6]  is  given  in  polar  coordinates  and  thus  requires  this 
additional  transformation.  Note  that  this  was  not  necessary  for  the  case  of  the  first 
expectation.  That  is  due  to  the  fact  that  the  polar  coordinate  radial  value  is  identical  to 
that  of  Yj^  since  Y^^  was  on  the  x-axis. 


Using  the  density  above  we  may  express  (6.82)  as 


= j j M(u,a)p(T2  <^i) 


T +00 


T +C30 


;r  cos 


= j |m {u,a) ^ ^exp 


0 r. 


S 

2-0 


a au 


-V  ■ y ■ sin 


2 2 2 
a -\-rQ  -cos 


2u  ■ cos'' 


V « y 


2/  y 


[a-yf 


cos 


(6.87) 


dadu 
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We  need  only  compute  the  final  expectation  now.  Just  as  with  the  first  two  terms, 
we  will  use  the  transformation  so  that  we  can  use  another  density  provided  in  [6],  Recall 
that  we  are  computing 


E9 


(6.88) 


We  break  this  up  into  two  separate  expectations 


e 


-rT 


+ e~^'^E9 


(6.89) 


We  focus  on  the  first  expectation.  First  we  remove  the  ( • Y notation  with  the  use  of  an 
appropriate  indicator  function.  This  yields 


-rT 


(6.90) 


Now  we  express  the  random  variable  in  the  expectation  in  terms  of  the  transformed 
process  F,  The  indicator  function  is  left  as  is  but  will  be  changed  momentarily. 

This  yields 


exp 


1 v2\  M\  tA 


r X 


) + c^i2 ))  ^i)j 


Again  we  simplify  the  notation  by  defining  the  function  by 


(6.91) 


.2  , ..2 


L{x,y)  = e "^exp  //,  ■ f [x,y)  + g{x,y)~  ^ Y x 

V 2 j (692) 

(^'  exp(o-„  •/(x,>')  + a,2 


Then  we  may  write  (6.91)  as 
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1.1 


(6.93) 


The  random  variable  in  this  expectation  depends  only  on  the  process  Yj.  when 

^3  = {r,  > T;  ^2  > r} . But  we  have  this  density  as  well  from  [6],  The  density  is  given  in 

polar  coordinates  as 


P(r,  > T,T2  > TJj  &dy) 


2r 

Ta 


exp 


f 2 , 2 ^ 

r + /•„ 


2T 


jn<j)  . jnG^ 


7=0 


sm  sm 
a 


a •'% 


V ^ y 


drd(f) 


(6.94) 


Just  as  before,  we  need  to  examine  the  indicator  function  a little  closer  in  order  to 
determine  the  region  over  which  we  will  integrate. 

{.sj.  -i,  2 0)  = {5'exp(cr,|  + so} 

( k ■) 

where  - In  ^ as  before.  If  we  write  (6.93)  as  a Lebesgue  integral  we  have 


(6.95) 


= j J leZ(rcos^,/-sin^)/’(z-,  >T,T2>  T,Y,  e drdO) 


(6.96) 


0 0 


where  •/(/'cos^,/-sin(9)  + CT,2  •g(rcos^,/-sin<9)  > Ajj.  We  may  rewrite  B as 
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B = {o-,j  (rcos^  • cos^  -r  sin^  -sin^  + Xg) 

+c7,2  {r cos ^ sin ^ + r sin (f> cos d +yo)>ky 
= {(/■  • cTji  cos G + r-  <Ti2  sin  cos <f> 

+ {r  ■ a,2  cos ^ - r • <7ji  sin sin ^ 

+r(c7„Xo+a,2>'o)^^i} 

= {('■•ctj2  cos^  - r c7,,  sin^)sin^  + 5 > 

= \r{(Ji2  cos^-cr,j  sin^)sin^  >^i  -a^ 

- sin  ^ - a| 


r > 


(^1  - <5) 


■ sin^ 

Define  the  function  (^)  by 


(6.97) 


^1  (<Z^)  = 


max 


k^-a 


7^ 


+ (Ti2  sin( 


(6.98) 


so  that  the  lower  bound  on  r for  a fixed  value  of  <t>  is  (^) . The  maximum  between  the 
quantity  above  and  zero  is  used  to  ensure  that  r > 0 in  the  case  when  k^<a  which  is 
equivalent  io  k^  <a . Now  we  may  express  (6.93)  as 


= 11  L{^rcos<f>,rsin</>)P[T^>T,T2>T,Y,Gdrdif>) 

0 9,(«l) 


(6.99) 


We  compute  the  second  expectation  in  just  the  same  manner.  We  write  out  the 
expectation  in  terms  of  the  transformed  process.  Recall  the  second  expectation  is 


[dP 

,9  , 

■€ 

19* 

1 

T ~ ^2) 

(6.100) 
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Written  out  in  terms  of  the  transformed  process  Y,  this  looks  like 


exp 


'4  Hfc,20|  ■ (■5'  «>‘P  K • / {rr . iV' ) + <^22  • « (n-‘. }?))  - *2 )) 

Define  the  fijnction  (/(■,■)  by 


(6.101) 


2 , ..2  A 


^/(x,>;)  = exp  //,  ■/(x,jy)  + ^.g(x,y)-^.'^^  T 


y 

(t?  exp (0-2,  ■f{x,y)  + a^  S {x, y)) - k^)"j 


(6.102) 


Then  (6. 101)  becomes 


(6.103) 


Again  the  random  variable  depends  only  on  the  process  Y,  when  = {r,  >T;  x^>  T] 
So  we  may  express  this  in  terms  of  the  Lebesgue  integral 


U ^ou 

= J J l^f/(rcos^zi,rsin^)2P(r,  >T,z^>  T,Y,  & dr  d<f>) 

00  ' 


(6.104) 


where  £ - {a,,  • /(/-cos<z>,/-sin^) + a,2  • g(/-cos^,r  sin^)  > ^’2}.  Calculations  similar 
those  of  (6.97)  show  that 


to 


E = 


\h-b) 


r \ 

yl^u  +^i2  J 


COS(p  + 


r \ 

^12^22 

Wn  + o-,2 


sin^ 


(6.105) 
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Define  the  function  q.^  (•)  by 


q2  {</>)  = max 


(*>*)■ 


+ 


(7,2<T2i  Q'iiO'22 

V^ll  ■*■^12 


0X2^22  ^11^21 


COS^ 


I 2 2 

■y/(Tii  +(Ti 


N 

-1  ^ 

sin^ 

,0 

y 

y 

(6.106) 


SO  that  the  lower  bound  on  r for  a fixed  value  of  (f>  is  ^2  (^)  • Then  we  may  express 
(6. 104)  as 


a .00  (6.107) 

= J I U [r  cos  (!>,r  sin  (!>)P{t^>T,  V2>T,Y,  & dr  d(f>) 

0 q2{<t>) 

This  being  the  last  necessary  calculation,  we  may  add  up  the  quantities  in  (6.73),  (6.87), 
(6.99)  and  (6  . 107)  in  order  to  obtain  the  price  of  the  get-out  option. 


in 


APPENDIX 

EXPLICIT  OPTIONS'  PRICING  FORMULAS 
The  results  for  the  partial  tunnel  options’  prices  from  chapter  4 were  expressed  i 
terms  of  the  functions  h(c,d),  h(c,d)  and  w(c,d)  from  chapter  3.  Here  we  present  the 

options’  prices  explicitly  in  terms  of  the  bivariate  normal  distribution  function.  We  first 
make  some  computations  that  will  be  used  later. 


Recall  that  = — 
a 


r 

V 2; 


So  the  quantity  ^ becomes 

2 


(/i  + cr)^r  +2o/i +CT^)r 
2 2 


+ 


f 

1 

f 

2 > 

2a 

r - 

a 

+ a^ 

V 

a 

\ 

2 J 

y 

T 


2 2 
_ ^ [2r-a^ 

" 2 ^ 2 
u^T 

= ^ — + /-T 


As  a result  of  theorem  4.2  we  have 

)-h{n  + a,dj))-k- [h{n,  ) - h{n,  dj ))]  (A.  1 ) 


-rT-^  n-T 

y=-cc 


where 
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h{c,d)  = exp 


-cd 


^ a-h  p-k^ 


-N^ 


^ a-d  p~k^ 


’ Vr 


yj 


Cj  - 2{b  - a) , dj  = 2{b  -a)-2b , p = , a=ct^-d  and  P = cT  -d . 

h{c,d)  and  simplifying  yields 

(p  + cr)/]  - Cj-d  (p  + a)T -Cj-k 


(A.2) 


Substituting 


+00 

S 

( 

-{ti+(T)c 

e ■' 

1 

N 

1 

V 

>/T 


-A, 


{fu  + a)t^  - c j - h {^  + a)T-Cj-k^ 


A 


Vf 


-e 


-(p+<T)i/, 


A. 


{^  + a)t^-dj -d  {^  + a)T -dj -k^'^ 


Vf 


y>' 


-A. 


(p  + <t)/i  - - Z>  {fi  + a)T-dj-k 


r\ 


7f 


Jy 


A. 


p/[  - - a fiT  -Cj  -ic^ 

7T  ’ y? 


-A. 


p/,  -Cj  -b  juT  -Cj  -k 

A ’ ^ . 


A, 


fit^-dj-d  i^T-dj-k^^ 


\ 


-N. 


yy 


^ /2t^-d  -b  /uT  -d  -ic 


jy 


. A ’ ^ 

Theorem  4.4  gave  us  the  price  of  the  PTPO-I  as 

K [k^, Cj )- Km, dj )) - *^ • (Km  + (T,Cj)-kM  + dj ))]  where 
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h{c,d)  = exp 


-cd 


^ b-a  k~P^ 


-N^ 


^ d-a  k - 


’ VF 


7J 


Cj  = 2(6  - a) , dj  = 2(6  -d)-2b , p = , a = ct^-d  and  P = cT-d 

Substituting  h(c,d)  and  simplifying  yields 


N 

'^b-pt,+Cj  ic-pt,+Cj^ 

^a-pT  + Cj  ic-pt,+Cj^ 

^ \F\  'Tt  ^ 

^p 

< 'Tt 

-S.{. 


-Hd 

9 J 


-iH+a)c^ 


N. 


-N_ 


N„ 


^ b-  pt^+dj  ic  - ///,  + dj  ^ 

Jr,  ’ ^/^ 

^ b - {p  + u)t,  + Cj  ic  ~{p  + a)T  + Cj^ 


''  a - pT  + dj  ic  - pt,  + dj  ^ 


^F^  ’ 'Tt 


yj 


-N„ 


^ d-{p  + cr)r,  +Cj  ic  ~{p  + a)T  + c.  ^ 


yj 


-{p+cj)d 

■>  J 


^ b ~{p  + a)t,+dj  ic  ~{p  + a)T +dj^ 


7F 


-iV. 


a-{p  + cr)t,  +dj  ic  ~{p  + a)T  + d- 


yj 


Theorem  4,5  gave  use  the  price  of  the  PTCO-II  as 


J^—<C 


-rT-U^T 


[5  • [w{p  + a,Cj)~  w{p  + o,  dj ))  - A:  • (w{p,  Cj  ) - w{p,  d^ ))] 


where 
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w{c,d)  = exp 


( c^T  ^ 
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r 

cd 

1 2 j 
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P 

\ 

b +d -a  b- P 
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^b+d-a  S-p^ 


-N^ 
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^d  + d-a  b~P^ 


n/^  ’ Vr 

^d+d~a  S~p 
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7^  ’ Vf 

5 = max(a,^),  = 2{b  - a),  dj  = 2(b-d)-2b , a = ct^-d  and  P = cT -d 

Substituting  w(c,d)  and  simplifying  yields 
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e ^ 
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b + 2c j - pl^  b - pT  + c. 

- N 
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P 

P 
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6 + Id.  - b - /iT  + d . 

. k ^ J 

^ b + Id^  - /i'/,  ^ - pT  + dj  ^ 
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-e  k ■ 
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7F 
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)) 


-N 


a + 2Cj  - b - /jT  + Cj^ 


- e 


-Md, 
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d + 2c^~  S - nT  + c. 
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b + 2dj  - b - fjT  + d 

b + 2dj  - /v/|  S - pT  + d 


^ ' Vf 


^ f d + 2d  - S - /jT  + d ^ 
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where  /u'  = p+cr  . 

Theorem  4.6  gave  use  the  price  Fjf  of  the  PTPO-II  as 


■ VF 


yj 


yy 


= z 


^ -rT-WT 


k ■ ( ) - w{p,  d^ ))  - .S  • (w^/u  +a,c^)-  w{^  + a,  d^ ))] 


where 
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w{c,d)  = exp 


r 

y 

y 
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1 2 J 
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b + d-a  y/ - P 
^ b + d ~a  d- 

~1T'^ 


-N. 
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•J<i  ■'VF 
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~K~’~W 


y/  = min(i,  k) , Cj  = 2(b  - d),  dj  = 2(b  - d)-2b  , a = ct^  - d and  = cT  - d . 

Substituting  w{c,d)  and  simplifying 

yields 
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